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1.0 INTRODUCTION 


1 . 1 PROGRAM GOALS 


This interim report describes computational methodology developed 
to analyze the plume flowfield generated by the interaction of an 
imperfectly expanded supersonic jet with the surrounding external 
stream. The overall goals of the present program are enumerated 
below. 


(1) the development of a "fully-coupled” parabolized Navier- 
Stokes (PNS) computer code to predict the multiple-cell 
wave structure in a two-dimensional (planar or axisymmetric) 
turbulent jet exhausting into a quiescent or supersonic 
external stream; 


(2) the formulation of techniques for the "strongly inter- 
active coupling" of the 2D jet solution with an external 
potential flow solution to provide for the analysis of 
jets exhausting into subsonic/transonic external streams; 


(3) the development of a fully-coupled PNS computer code to 
analyze the nearfield structure of three-dimensional 
turbulent jets (i.e., jets issuing from rectangular 
nozzles) exhausting into a quiescent or supersonic ex- 
ternal stream; and, 


(4) the formulation of techniques for the strongly interactive 
coupling of the 3D jet solution with an external 3D 
potential flow solution. 


This report discusses the accomplishments achieved in satisying 
the first goal listed above. Part II of this report describing 
accomplishments achieved in satisfying the second goal is now in 
preparation. The 3D methodology of the last two goals will be a 
direct extension of the 2D methodology. This 3D work is now in 
progress and will be documented in a forthcoming report. 


1.2 APPLICATIONS 


This effort is jointly supported by the Propulsion Aerodynamics and 
Aeroacoustics Branches of the NASA Langley Research Center , whose 
applications for this technology are quite disparate. The Propulsion 
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Branch is concerned with the influence of the jet exhaust on the 
nozzle afterbody pressure distribution for subsonic/transonic 
flight conditions. In a previous effort, Dash and coworkers 
developed an overlaid viscous/inviscid jet model 1 " 2 * 4 which was 
integrated by Wilmoth 5 * into a patched component nozzle afterbody 
model. This methodology was demonstrated to work quite well in 
weakly interactive situations 6-8 , but, could not reliably treat 
strong viscous/inviscid interactions, and, was not readily ex- 
tendable to three-dimensional jet flowfield problems. In this 
program, new fully-coupled technology to deal with strongly in- 
teractive jet phenomena has been developed 9 * 1 0 which is readily 
extendable to 3D flow problems, as now being performed. 


The Aeroacoustics Branch is concerned with the prediction of shock 
noise in imperfectly expanded jets. This requires a detailed 
portayal of the coupled multiple-cell wave/shock and turbulent 
mixing processes occuring in the jet (see references 11 - 13) which 
has not heretofore been available. The technology formulated in 
this program has led to the development of a new model and asso- 
ciated computer code 14 which has been exhibited 15 to provide this 
capability . 


1.3 ACCOMPLISHMENTS 


The accomplishments achieved to date in this program are summarized 
below. 


(1) A 2D (planar/axisymmetric) PNS jet model (SCIPVIS) 14 has 

been developed which calculates the fully-coupled viscous/ 
inviscid jet interaction flowfield for a jet exhausting 
into a quiescent or supersonic external stream. A single- 
pass explicit spatial marching procedure is used to 
perform this calculation. SCIPVIS combines hyperbolic/ 
parabolic shock-capturing methodology for treating super- 
sonic portions of the jet with partially-parabolic 
methodology for treating subsonic portions. Compressibility 
corrected two-equation turbulence models are utilized to 
determine the local turbulent diffusivity, and specialized 
techniques are employed to locate Mach discs and calculate 
the interactive subsonic/super sonic mixing region downstream 
of the Mach discs. 


(2) A series of ducted supersonic mixing calculations were 

performed to exhibit the ability of SCIPVIS to analyze 

fundamental interactive phenomena 9 * 10 (i.e., waves gen- 

erated by supersonic mixing and wave/shear layer inter- 

actions ) . 
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(3) A series of under and overexpanded free jet calculations 

were performed for jets exhausting into both quiescent and 
supersonic external streams to exhibit the performance of 
the SCIPVIS code under a broad range of operating 
conditions. The quiescent stream calculations were com- 
pared with the laboratory data of Seiner and Norum 11 in a 
code assessment effort performed by Seiner at NASA/LRC, 
which is documented in reference 15. 


(4) A 2D partially-parabolic implicit model (SPLITP) 9 has been 
developed which calculates subsonic/transonic wall bounded 
and free mixing regions. The wall bounded version of 
SPLITP 16 * 17 was developed for David Taylor Naval Ship R&D 
Center to analyze curved wall jets over the aft portion of 
circulation control airfoils. For jet/afterbody inter- 
action applications, SPLITP is applied in the wall bounded 
mode over the nozzle afterbody and in the free jet mode, 
downstream. A direct coupling approach has been formulated 
for coupling SPLITP with an external potential flow solver 
which utilizes pressure-splitting methodology in non-reversed 
flow regions. A velocity-split coupling approach has been 
proposed for analyzing reversed flow regions. SPLITP pres- 
ently serves as a starid-along code for use in investigating 
jet/potential flow coupling techniques. A description of 
SPLITP and the coupling methodology developed will be pro- 
vided in Part II of this report. 


An envisioned final product of this technology for j et/afterboay 
applications would be a unified get/boundary layer model which 
combines the SCIPVIS and SPLITP models. The SPLITP portion would 
analyze the afterbody boundary layer and the subsonic portion of 
the jet mixing layer ; the SCIPVIS portion would analyze the 
supersonic portion of the jet mixing layer. Direct pressure-split 
coupling would be utilized in moderately interactive situations 
while velocity-split coupling would be employed in strongly 
interactive situations . This envisioned approach is conceptually 
the same for 2 or 3 dimensional problems. 


This report (Part I) details the computational features of the 
SCIPVIS model and present results obtained for a spectrum of free 
and ducted problems. The SCIPVIS and SPLITP computer codes have 
been delivered to NASA/LRC and are operational at that facility. 
Separate documentation describes the operation of these codes. 
Part II of this report (Pressure-Split Model, SPLITP) describing 
the computational features of the SPLITP model and results for a 
number of free jet and wall jet problems is now in preparation. 
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2.0 OVERVIEW OF FULLY-COUPLED METHODOLOGY 


2.1 JET FLOWFIELD CHARACTERISTICS 


The analysis of the plume flowfield generated by the interaction of 
an imperfectly expanded jet with the surrounding external stream 
involves consideration of a number of distinct flow regions with 
varying characteristic features and length scales. Figure 1 exhibits 
the nearfield features of an under expanded jet exhausting into: 


(1) a quiescent external stream; 

(2) a supersonic external stream; and, 

(3) a subsonic external stream with a large boundary layer 
present . 


In all these situations, the jet portion of the nearfield is charac- 
terized by a predominantly inviscid shock cell structure with mixing 
layers growing along the jet and Mach disc slipstreams. A transi- 
tional region (Figure 2) joins the predominantly inviscid nearfield 
with the fully viscous pressure-equilibrated farfield. Here, the 
mixing layers come to engulf the entire jet and wave processes occur 
in a fully turbulent environment. Wave processes are confined to 
supersonic regions of the flow bounded by the viscous sonic lines in 
the jet and/or Mach disc mixing layers (see Figure 1). 


In the jet nearfield, strongly interactive phenomena occur in many 
situations which cannot adequately be analyzed using the earlier 
patched component (overlaid) methodology 1- *'. Such interactive 
phenomena include: 


(1) compression waves produced by the positive displacement 
effect of chemistry or high Mach number viscous dissipation 
in the jet mixing layer; 

(2) expansion waves generated by the "washing away" of large 
mass defect initial regions downstream of base/separated 
flow zones; 
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{• FIRST SHOCK CELL 


• QUIESCENT EXTERNAL STREAM 



• SUPERSONIC EXTERNAL STREAM 



• SUBSONIC EXTERNAL STREAM 


FIGURE 1 

Schematic of Jet Nearfield Structure for Quiescent, 
Supersonic and Subsonic External Flows. 



• DAMPING OF WAVE AMPLITUDES 



FIGURE 2 

Schematic of Jet Nearfield, Farfield and Transitional 
Region; Damping of Wave Amplitudes by Turbulence in 
Transitional Region. 


(3) the negative displacement of streamlines downstream of 
Mach discs (see Figure 1) generated by the wave-like 
entrainment process occurring; and, 


(4) shock/shear layer interactions occurring at the end of 
shock cells for situations where a significant portion 
of the nearfield jet core has been entrained into the 
jet layer. 

The above interactive phenomena are schematized in Figure 3. 


The transitional region (Figure 2) where wave processes are large 
ly embedded in the turbulent mixing layer, is always a region of 
strong viscous/inviscid interactions and cannot be treated by 
weakly interactive patched component technology. Here, wave 
intensities are damped by turbulent dissipation and wave fronts 
are curved by the higher mixing layer rotationality as exhibited 
in Figure 2. For many engineering applications (i.e., plume sig- 
nature predictions or plume/af terbody interaction problems), the 
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• SOURCE -LIKE DISPLACEMENT - PRODUCED BY TURBULENT DISSIPATION 



• SINK-LIKE DISPLACEMENT -PRODUCED BY MIXING DOWNSTREAM OF 




• WAVE/SHEAR LAYER INTERACTION - TRANSMISSION, REFLECTION, 



Fundamental Interactive Phenomena Occurring in Supersonic 
Mixing Flowfields. 
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decaying wave structure in the transitional region has a negligible 
influence and can be approximated. For the analysis of jet 
laboratory data or for applications to problem areas such as jet 
shock noise 11 *" 13 , the details of the shock structure in the trans- 
itional region are quite important. 


The above interactive phenomena are all "quasi-parabolic" in the 
sense that diffusive processes along streamlines can be neglected. 
This report discusses several quasi-parabolic techniques developed 
to analyze such phenomena via the solution of a fully-coupled 
parabolized system of viscous/ inviscid equations and the use of 
"direct-coupling" techniques to join together the various regions 
of the plume. The analysis of plume induced separated flow regions 
requires the use of fully-elliptic viscous/inviscid methodology 
which has not been formulated in this effort. An extension of the 
quasi-parabolic methodology using velocity-split coupling techniques 
can perform this type of analysis in significantly less time than 
full Navier-Stokes methodology and is now under investigation at 
NASA/LRC. An overview of the extensions required will be provided 
in Part II of this report. 


2.2 BASIC MODELING REQUIREMENTS 


The prediction of interactive jet flowfields entails the unification 
of modeling techniques for analyzing the following processes: 


(1) wave/shock propogation in both inviscid and viscous 
supersonic regions of the jet; 


(2) turbulent mixing in the nearfield jet shear layer 

(Figure 1) and in the transitional region (Figure 2); 


(3) the influence of compressibility effects and pressure 
gradients on the turbulence; 


(4) the strong interaction of the wave and turbulent 
structure as occur at the end of shock cells and 
throughout the transitional region; 

(5) the interactive coupling of the viscous/inviscid jet 
with a subsonic or supersonic external stream; 
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(6) the occurrence of Mach discs and the strong interactions 
induced by the wake-like turbulent mixing process occurr- 
ing behind the disc; and, 


(7) subsonic/supersonic coupling at the viscous sonic lines 

occurring in the jet and Mach disc mixing layers (Figure 1). 


2.3 ANALYSIS OF WAVE/SHOCK PROCESSES IN SUPERSONIC INVISCID AND 
VISCOUS REGIONS 


The treatment of wave/shock processes in inviscid flow regions is 
performed utilizing the SCIPPY model shock-capturing approach of 
Dash and Thorpe 1 e “* 9 . Here, the equations are cast in conservative 
form and integrated utilizing the alternating one-sided difference 
algorithm of MacCormack 20 which permits embedded shock waves to be 
"numerically captured". This approach has yielded results equiva- 
lent to those obtained utilizing more cumbersome shock-fitting 
algorithms for a variety of jet exhaust flowfield problems (see 
Reference 19). 


The extension of this shock-capturing approach to supersonic viscous 
flow regions involves the addition of parabolized stress/transport 
terms which render the resulting viscous/inviscid jet equations 
hyperbolic/parabolic. The use of a fully-coupled system of viscous/ 
inviscid equations to analyze wave processes in mixing regions was 
first introduced by Ferri in the mid-sixties for application to 
supersonic combustion problems (see the review article by Ferri 21 for 
details of this earlier work). The numerical solution of these 
original fully-coupled equations employed a "viscous-characteristic" 
procedure 22 and a wide spectrum of problem areas were treated 
utilizing this technology. In present applications, the fully- 
coupled equations are termed the parabolized Navier-Stokes (PNS) 
equations and a number of numerical algorithms have been developed 
for their solution (see the review article of Dash and Wolf 9 ). The 
shock-capturing algorithm utilized in the SCIPVIS model is described 
in detail in this paper. The ability of this shock-capturing algo- 
rithm to analyze basic interactive phenomena occurring in supersonic 
mixing regions 10 (i.e., the interaction of expansion fans and shock 
waves with shear layers, waves generated by high Mach turbulent 
dissipative processes, etc.) is demonstrated by numerical studies of 
ducted mixing. Complete jet flowfield solutions for supersonic jets 
exhausting into supersonic streams, predicted using this shock- 
capturing methodology, will be discussed. These studies have delin- 
eated conditions under which weakly interactive overlaid viscous/ 
inviscid coupling 3 becomes inapplicable and fully-coupled methodology 
is required. 
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2.4 ANALYSIS OF QUASI -PARABOLIC SUBSONIC VISCOUS REGIONS 
2.4.1 Embedded and Adjoining Subsonic Regions 


In the interactive supersonic jet problem, we encounter both embed- 
ded and adjoining regions of subsonic quasi-parabolic viscous flow. 
Defining the jet portion of the flowfield to be that portion of the 
overall interactive flowfield bounded by the jet axis and the outer 
edge of the jet mixing layer (see Figure 1), two types of embedded 
subsonic regions occur, namely: 


(1) the region behind Mach discs (Figure 1); and, 

(2) the initial boundary layer dominated region of the jet 
mixing layer for jets exhausting into a supersonic 
stream (see Figure 2 - sinklike displacement schematic). 


Adjoining subsonic viscous regions occur for supersonic jets ex- 
hausting into quiescent or subsonic/transonic external streams and 
occupy the region of the jet bounded by the jet mixing layer sonic 
line and the mixing layer outer edge (see Figure 1). 


2.4.2 Pressure-Splitting Approximation 


The governing flowfield equations in these subsonic regions are the 
same fully-coupled PNS equations used to analyze supersonic viscous 
flow regions. However, due to the elliptic character of the flow, 
a different numerical approach is required to integrate these 
equations. The approach taken here permits spatial marching in sub- 
sonic regions via the splitting of the pressure field such that the 
streamwise pressure gradient is "stipulated" while the crossflow 
pressure variation is arrived at via the coupled solution of the 
continuity and normal momentum equations. This type of approach has 
been implemented for both ducted and free jet mixing problems by 
Spalding and coworkers 2 3 » 2 4 , Briley and McDonald 2 5 » 2 6 » and a number 
of other investigators. 


Different types of splitting approximations are used to deal with 
the different subsonic flow regions. Downstream of Mach discs, the 
streamwise pressure gradient is suppressed in a manner akin to that 
implemented in supersonic external flow PNS solutions in the near- 
wall region 27 (i.e. , a "sublayer" approximation is utilized). 
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This approximation presumes that the Mach disc radius is small and 
hence, the acceleration of the subsonic flow behind the disc to 
supersonic velocities is dominated by the turbulent mixing process. 
The transverse pressure variation across the Mach disc mixing 
region is neglected under these circumstances. In analyzing the 
subsonic portion of the jet mixing layer for a quiescent external 
stream, the streamwise pressure gradient is set to zero; for sub- 
sonic/transonic external streams, it is set equal to the pressure 
gradient existing at the outer mixing layer edge as established by 
a potential flow solution. 


2.4.3 Direct -Coupling of Jet and Potential Flow Solutions 


The jet and potential flow solutions are "directly-coupled" at the 
outer edge of the jet mixing layer utilizing an approach that para- 
llels that of Bradshaw and coworkers 2 8 > 2 8 developed for boundary 
layer /potential flow coupling. Implementation of this approach has 
involved : 


(1) solution of the jet equations in a coordinate system 
which maps the outer edge of the jet mixing layer to 
a constant in mapped coordinates; 

(2) prediction of the jet boundary variation via an ordinary 
differential equation which accounts for both jet en- 
trainment and plume expansion; 

(3) solution of the subsonic jet equations using pressure- 
split methodology with the streamwise pressure gradient 
stipulated along the jet outer boundary; 


(4) revision of the streamwise pressure gradient via solution 
of the potential flow equations over a geometric surface 
comprised of the jet boundary with a transpiration 
boundary condition corresponding to the entrainment 
(radial) velocity component at that position; and, 

(5) the use of a "global" pressure update procedure 28 ’ 28 
involving several sweeps of Step (3) for a fixed outer 
streamwise gradient to analyze regions with significant 
streamwise and/or normal pressure variations. 
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Steps (1) - (3) are operational in both the SCIPVIS and SPLITP 
computer codes. Steps (4) and (5) are under development using 
SPLITP. The details entailed in performing the above operations 
will be provided in Part II of this report. 


2.4.4 Subsonic/Supersonic Coupling 


In analyzing mixed subsonic/supersonic viscous flows, a formal 
coupling procedure is required to join the two solutions smoothly 
at the viscous sonic lines in the jet and Mach disc mixing layers 
(see Figure 1). The approach employed parallels that developed by 
Ferri and Dash 22 to couple the subsonic and supersonic portions of 
an interactive supersonic boundary layer solution. In this approach: 


(1) convective and diffusive derivatives at the matching 
grid point smoothly connect variables on either side; 
and, 


(2) a "viscous-characteristic" compatibility relation on 

the supersonic flow side ensures compatibility between 
the pressure and flow angle at the matching point. 


In situations where normal pressure variations across the subsonic 
viscous layer are analyzed, and/or where the streamwise pressure 
variation is locally determined via an integrated continuity con- 
straint (as in analyzing the Mach disc mixing region), this coupling 
procedure is iterative. 


2.5 TURBULENCE MODELING CONCEPTS 


Turbulent mixing processes are represented using classical Boussinesq 
type approximations bo relate the turbulent shear stress and scalar 
transport terms to the mean flow gradients. The turbulent diffusivity 
is determined using conventional models of the two-equation class 
which solve partial differential equations for the variation of the 
turbulent kinetic energy and a length scale parameter. This level 
of turbulence modeling represents the present state-of-the-art for 
"practical" analyses of high Reynolds number jet f lowf ields 3 0 . Ap- 
plications to a complete spectrum of simple (constant pressure) free 
jet problems 31 had indicated that basic two-equation turbulence 
models performed reliably for low speed jet problems but performed 
poorly in supersonic mixing situations. Present generation two- 
equation models are all formulated from an incompressible viewpoint 
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and thus, their inability to deal with such compressible phenomena 
is not surprising. Extensions of two-equation models to deal with 
compressibility effects have thus far been heuristically formulated 
(i.e., no attempt has yet been made to formally model the terms 
embodying the compressibility effects of high Mach numbers). A 
detailed assessment of the performance of compressibility modified 
two-equation turbulence models for both simple and complex (variable 
pressure, chemically reacting) jet flowfields is given in reference 
30. An overview of high Mach number compressibility effects on 
turbulent mixing and the specific "compressibility corrected" tur- 
bulence models employed in this effort will be discussed in this 
report. No attempt has been made to model the large scale structure 
of the turbulence 32 associated with the acoustic excitation of jet 
instability modes. 
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3.0 GOVERNING EQUATIONS 


3.1 REYNOLDS DECOMPOSED JET MIXING EQUATIONS 


The time-averaged, Reynolds decomposed, planar (J=0) or axi- 
symmetric (J=l) parabolized jet mixing equations 33 are listed 
below. 
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These equations have been parabolized with respect to the axial 
(x) direction (viz. all transport terms with axial derivatives 
have been deleted) and "standard" (incompressible) assumptions 
concerning the turbulence correlations have been invoked (v iz . , 
thir d and hig her-order correlations, density correlations {p'u', 
p ' H ' and p ' 4> ' } etc. are neglected). The transport of heat and 
mass is taken to be the same (i.e., the Lewis number is taken to 
be unity) and thus only the laminar/turbulent Prandtl numbers 
appear in the energy and species transport equations. In the 
above equations, U and V are the axial and radial velocity com- 
ponents, p is the density, P is the pressure, p is the laminar 
viscosity, H is the total enthalpy, and <(> is the species para- 
meter. 


3 . 2 THERMODYNAMICS 


The jet mixing problems considered assume that the jet and external 
streams are each of uniform composition. For nonreacting (chem- 
ically frozen) situations, the species parameter 4> describes the 
local mixture composition; viz., 



( 6 ) 


where ai is the mass fraction of the ith species and J and E rep- 
resent the constant values of ai in the unmixed jet and external 
streams. The static enthalpy is given by: 


h( 4> , T ) = (h J (T) - h E (T)}<{> + h £ (T) 


where : 


J , E 


= E {a h (T) } 


i i ' J,E 


Then the specific heat ratio, y($.T) is given by: 


( 7 ) 


C p (<|>,T)W(4>)/R 0 
Y( * ,T) “ Cp( <j>,T)W(<|))/R 0 -l 
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( 9 ) 


where the specific heat, 



is given by: 


C p (*,T) 


3h 

3T 


and the molecular weight, W, is given by: 
3.3 PARABOLIZED STRESS TERMS 


Eliminating all terms containing streamwise derivatives, the 

laminar stress terms x , x and x are as follows : 

xr xx rr 


T xr 

= up 

w 3r 

(Ha) 

T xx 

= _ 2 JV 
3 y 3r 

(Hb) 

T rr 

= ±J1 

3 y 3r 

(He) 


Using the Boussinesq-type approximation: 


- pu . 'u ' . 
i J 



c 


3Ui 
3X . 





( 12 ) 


(where the turbulent kinetic energy k “m'ui'/ 2 )* the para- 
bolized turbulent stress terms are as follows: 


pu'v' 

-Uti? 


(13a) 

pu'u' 

- - l pk 

2 3V 
“ 3 y t3r 

(13b) 

pv' v' 

= - i pk 

4 3V 
3 y t3r 

(13c) 
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The turbulent transport of scalar variables, ot, is expressed by.: 


- pa ' v ' 


^t 9a 
a 9r 


a 


( 14 ) 


where a represents the turbulent Prandtl number, Pr t , for a = H 
or <f>. a 


3.4 TWO EQUATION TURBULENCE MODELS 


In the generalized Prandtl-Kolmogorof f turbulence formulation 
(see, e.g., reference 34) the turbulent viscosity, Pt> * s related 
to the turbulent kinetic energy, k, and turbulent length scale, 
l, via: 


U 


t 


c y pA 


(15) 


where Cy is a dimensionless constant. A partial differential 
equation can be formally derived to describe the production, dis- 
sipation, and transport of turbulent kinetic energy. By analogy 34 , 
an equation can be formulated for the length scale parameter, z, 
which is related to k and & by: 

z = k m £ n (16) 

Table 1 shows three different popular forms of the length scale 
equation . 


Table 1 - Forms of the Length Scale Equation 


Model 

z 

m 

n 

W 

ke (Launder, et al.) a 

E 

3/2 

-1 

k 2 /e 

kW ( Spalding ) b 

W 

1 

-2 

k/W* 

km (Saffman) 0 

10 

1/2 

-1 

k /to 


a See Reference 31 
^ See Reference 35 
c See Reference 36 
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Only the ke 31 and kW 35 two equation models have been incorporated 
in SCIPVIS and SPLITP. The kia model of Saffman 36 (in particular, 
the modified version which accounts for compressibility 37 ) does 
show promise, as demonstrated by Walker 38 in his evaluation against 
Mach 2.2 jet data. However, this model has not been assessed 
against a singificant body of compressible jet data. 


3.4.1 ke Turbulence Models 


The ke models 31 incorporated in SCIPVIS and SPLITP (kel and ke2) 
are extended forms of the basic model which contain both axisym- 
metric and weak shear flow "corrections." The following equations 
are solved for k and e : 


jZ (pUkr J ) + JL (pvkr J ) 
JL ( pUer J ) + (pVer J ) 



(17) 


(18) 


+ rJ ( C lE - C 2 pe) i 


In the above equations, P is the turbulent production term: 


- y t^8r^ 


(19) 


The turbulent viscosity, ut * is determined from the local values 
of k and e via the relation 

k 2 

V t = CjjCf.gJpf (20) 

The following model constants 31 are utilized: 



kel 

ke2 

°k 

1.0 

1.0 

°e 

1.3 

1.3 

C 1 

1.43 

1.40 

C 2 

1.92 - . 0667 f 

1.94 - . 1336f 

C 

u 

.09 - . 04f 

. 09g(P/e ) - . 0534f 
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The 

the 

and 


axisymmetric correction parameter, f, is set equal to zero in 
shear layer region (before the mixing zone reaches the axis), 
is defined by 



( 21 ) 


downstream of that point (r„ is the width of the full mixing layer, 
U c is the jet centerline velocity and U e is the external stream 
velocity). The weak shear flow correction, g, is a function of the 
shear-stress weighted average ratio of production to dissipation 
rates. The functional form of g(P/e) is given in Figure 4. 



P/e 


FIGURE 4 


Variation of g With P/e For ke2 Model 
(from Launder, et al . l ) 
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3.4.2 Compressibility Corrected Version of ke Models (ke,cc) 




The performance of the ke models in analyzing supersonic jet and 
shear layer data has been shown to be quite poor^ 1 . An effort 
was undertaken to correct this model in a heuristic fashion 39 to 
account for the reduced mixing rates observed for higher Mach 
number jet mixing. The compressibility corrected viscosity is 
given by: 


y t = K < W T 


( 22 ) 


where K(M T ) is the correction factor and M T is the characteristic 
Mach number of the turbulence (M x = /kmax/a, where k ma x is the 
maximum value of k at each station and a is the local sound speed 
at the grid point where k is maximum). The functional form of 
K(M X ), given in Figure 5, was determined by matching calculations 
to observed spread rates for isoenergetic , supersonic shear 
layers with one stream stationary 40 . Figure 6 exhibits the per- 
formance of the ke2 and ke2,cc models in predicting this spread 
rate data; the ke2 model shows no variation with Mach number while 
the ke2,cc model duplicates the data (as per the calibration of 
K(M X ) ) . Note that the compressibility correction term goes to 
unity as compressibility effects diminish and hence, the ke2 and 
ke2,cc models are equivalent in lower speed situations. 


3.4.3 kW Turbulence Model 


The kW model 35 , although developed at about the same time as the 
ke model, has not been widely used due to complications in applying 
it to wall-bounded shear flows 41 . However, a new version of the kW 
model 42 with coefficients set by rocket plume data has been 
employed successfully by the Propellant and Rocket Motor Establish- 
ment (PERME) in Great Britain. Predictions by PERME simulating 
laboratory rocket plume experiments performed at AEDC 4 3 were quite 
promising and the overall assessment of this model by Pergament 30 
for a spectrum of jet mixing problems was favorable. Hence this 
model has been incorporated into SCIPVIS and SPLITP. 
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Spreading parameter, 



The following equations are solved for k and W: 


(pUkr J ) + fp (pVkr J ) = ^ 


9r 


( 


r y t 3k. 
a, 3r ' 


+ r J (P - pC d kW^) 


fe (fiOTrJ ) + 1? <f> VKrJ > - 1? < 


+ 4 


C 1 W - “ C 2 pkW 


3/2 


r y 


'W 


t 3W. 
3r' 


+ C 3 y 


t 

1 3r^ J 


( 23 ) 


(24) 


The turbulent viscosity is given by 


y 


t 



(25) 


and the following model constants 42 are utilized 


a k = * 86 

a = .86 

e 

C D = .09 

C ± = 1.48 
C 2 = .18 

C 3 = 3.5 


Note that this model contains no axisymmetric or weak shear flow 
correction terms. 

3.5 SCIPVIS MAPPED, VECTORIZED CONSERVATION EQUATIONS 


The system of mean flow and turbulence model equations are solved 
in weak conservation form in SCIPVIS utilizing the mapped coordinates 
E, and n given by the simple rectangular transformation. 
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€ » x 

n = (r-r L (x))/(r u (x)-r L (x)) 


( 26 ) 


where tl and ry are the boundaries of the jet domain being solved. 
With this transformation, the equations can be expressed in the 
following vector form: 


3E 3F b^ J_ | r y eff 3f \ F 

3C an r J 3n l a f 3n / f 


(27) 


where : 


- [l, U, V, H, *, 


k (1 >, e , k< 2 > 


(1) designates ke model 

(2) designates kW model 




E = 


e p 


PU 


PV 

e u 


aP + pU 2 


pu\ f 

e v 


PUV 


P + pV 2 

e H 


pUH 


pVH 

e $ 

e (1) 
e k 

** 

pU4> 

pUk (1) 

; F = b 

PV4> 

pVk (1) 

e e 


pUe 


pVe 

V 2) 


pUk (2) 


pVk (2) 

_ 


pUW 


pvw 


-aE 
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and 


where : 



0 

(ct-1) 9P/ 9x 
0 

b 2 9 / r *^ef f ( a H -1 ) _9_ ,U^\ 

r J 3r V °H ^ 2 7 

0 

P - e 

(C^P. - C 2 pe) e/k 


P - pC D kW^ 

3 

CjWP - C 2 pkW T + C 3 P t 
- 

In the above equations, a and b are the transformation parameters 
given by : 



a(£,n) = (l-n)r^ + nr^/Cry - r L ) 
b(C) = l/(r D - r L ) 


and a is a flow indicator (a = 1 for supersonic flow regions, a = 

0 for subsonic flow regions). The variable V e ff is the combined 
laminar and turbulent viscosity. In practive, Ueff is set equal 
to the turbulent viscosity for high Reynolds number jets and the 
turbulent stress terms in the radial momentum equation are deleted. 
The analysis of lower Reynolds number jets would entail the use of 
a low Reynolds number turbulent model formulation. Low R e two- 
equation turbulence models have been developed for boundary layer 
problems to treat laminarization and the near wall region. Their 
ability to predict laminar to turbulent transition is questionable 
(generally, an empirical transition relation is incorporated and a 
spot of turbulence is artifically introduced) and the coefficients 
have not been calibrated for free jet mixing problems. Hence, 
while low Reynolds number jet problems can be analyzed, substantial 
work in extending present turbulence models formulations may be 
required . 
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3.6 VISCOUS CHARACTERISTIC EQUATIONS 


Boundary points in supersonic flow regions and subsonic/supersonic 
matching points are analyzed in SCIPVIS using formal characteristic 
relations. In viscous regions, the transport terms are included 
in the characteristic relations as local source terms using a 
viscous-characteristic formulation 22 . Along a X ± characteristic 
(Mach wave) given by: 


X 1 = = tan( 6 ± y) 


(29) 


the following pressure/flow deflection angle ( 0 ) compatibility 
relation applies: 


. giniI _ cp8M + d0 
Y 




-JsinB 

r 


f„ *i . ~ 

V siny 

^ 7 - dx 

yPM cos (0 ± y) 


(30) 


In the above relations, y is the Mach angle (siny 
is the viscous source term given by: 


where; 


F V = S Q - (Y-l) 


» 2 [? - ^ 


o _ cose d_ ( j 
S Q " 9r l r y eff 


'r J y 

Q _ cose S_ I M ef f 
H J 9r 


(* J If) 

( r y eff iH\ 

V o H 9r/ 




S_ f r<Jy eff v 
9r 


H 

- 1 ) 


H 


1_ (Hi) 

Sr '‘2 ; 


2 2 i 

and Q is the total flow velocity (Q = (U +V ) 2 ) 


1/M) and F^ T 


(31) 


C32a) 


C32b) 
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4.0 NUMERICAL PROCEDURES IN SCIPVIS 


4.1 INTEGRATION OPTIONS 


SCIPVIS provides for the spatial integration of the conservative 
jet mixing equations described by equation (27) using four types 
of integration options whose features and areas of applicability 
are described below. 


4.1.1 Parabolic (Constant Pressure Mixing) Option 


This option can be utilized to analyze balanced pressure subsonic 
or supersonic jet mixing. In utilizing this option: the pressure 

parameter a in equation 27 is set equal to zero (which removes the 
pressure term from the streamwise momentum conservation variable, 
eu)‘, the normal momentum equation is eliminated; the streamwise 
pressure gradient term, 9P/8x, in the Gf vector array is set equal 
to zero; and, a parabolic decode procedure (to be described below) 
i.s used to extract the mean flow? variables, U, H, and <f> from the 
E conservation variables, eu, eg, a nd In analyzing a balanced 

pressure jet (Figure 7), the equations are mapped from the jet 
axis to the outer edge of the jet mixing layer, rp(x ) , along which 
the external flow conditions Uj, Hj, and <|> = 0 are prescribed. 



SCIPVIS PARABOLIC (TURBULENT MIXING) MODE 


FIGURE 7 

Parabolic (Constant Pressure Mixing) Option of SCIPVIS for 
Analyzing Balanced Pressure Turbulent Jet Mixing. 


26 



The inner mixing layer boundary "floats" across the grid points 
until it reaches the axis. The parabolic boundary growth rate, 
drp/dx, is based on the outer edge gradients of streamwise 
velocity, U, and species parameter, <f>. and is given by: 


ffSl - (il) 

dx VIS f L 1 8r l 


( 33 ) 


where f represents U and 4> (the maximum gradient is utilized) and 
a value of C ' 1 places the computational boundary, rp» in close 
proximity with the physical mixing layer boundary. Increasing C 
yields a buffer region of nonturbulent flow between the physical 
mixing layer edge and rjj but does not effect the calculation except 
for decreasing the grid resolution. This integration option is 
employed in the analysis of underexpanded jets exhausting into a 
quiescent stream for analyzing the subsonic portion of the jet 
mixing layer, as will be described below. 


4.1.2 Partially Parabolic (Pressure-Split) Option 


This is an extension of the above option to account for streamwise 
and/or normal pressure variations in. the subsonic portion of the 
jet mixing layer using a pressure-splitting approximation. Here, 
the above parabolic option is first exercised (in integrating the 
equati ons fr om x to x + Ax) with the streamwise pressure gradient 
term, 8P/9x, imposed*. This yields the mean flow variables U, H, 
and <f> at x + Ax. 


The pressure, density, and radial velocity variation across the 
mixing layer are then arrived at via the coupled integration of 
the continuity and normal momentum equations and the equation of 
state constraint listed below (for a perfect gas); 


* The procedure for determining the imposed streamwise pressure 
gradient differs in accordance with the subsonic region being 
analyzed, as will be discussed. A detailed review of pressure- 
split methodology will be given in Part II of this report. 
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H 


(34) 


_1_E + i (u2 + v *> 

where H and U are known from the parabolic integration and thus, 

P, p, and V at each grid point are related in a nonlinear fashion. 

In the standard (constant pressure) parabolic mode, the continuity 
and above state relation are solved in a coupled fashion to yield 
the variation of radial velocity and density across the jet mixing 
layer. The details of the variable pressure crossflow integration 
will be given in Part II of this report which describes the numerical 
procedures in SPLITP. This option has been exercised and employed 
in SPLITP but has not been implemented in the SCIPVIS studies per- 
formed to date. The partially parabolic option is also implemented 
to analyze the subsonic turbulent mixing region behind Mach discs and 
will be described in a later section of this report. 


4.1.3 Hyperbolic (Inviscid Supersonic) Option 


This option can be utilized to predict the nearfield inviscid wave/ 
shock structure in the jet (as a standalone option, it reduces to 
the inviscid SCIPPY model 18 * 19 ), and is, of course, utilized to 
predict the inviscid core region of underexpanded jets. To predict 
the inviscid supersonic jet structure (Figure 8): the transport 

terms in equation (27) are set equal to zero; the turbulence, equa- 
tions are eliminated; and, the equations are integrated with the 
parameter a set equal to one (i.e., the pressure is included in the 
eu conservative variable). The inviscid jet equations are mapped 
from the axis to the jet interface, ru, defined by: 



(35) 


In performing an inviscid jet calculation, the pressure along the 
interface, r U( would be prescribed for quiescent or subsonic 
external flow and determined by a concurrent external flow solution 
(or the use of pressure/flow deflection rules such as shock expan- 
sion theory) for supersonic external flows (see reference 19 for 
details). A supersonic decode procedure is used to extract the 
flow variables: U, V, P, H, <j>, and p from the conservation vector, 

E. 
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/ 

SCIPVIS HYPERBOLIC (SUPERSONIC WAVE) MODE 


FIGURE 8 

Hyperbolic (Inviscid Supersonic) Option of SCIPVIS for 
Analyzing Inviscid Shock/Wave Structure in Underexpanded 
Jet . 


4.1.4 Hyperbolic/Parabolic (Viscous Supersonic) Option 


This option is utilized to analyze supersonic viscous (turbulent) 
regions of the jet. It is a direct extension of the hyperbolic 
option with the transport terms evaluated and the turbulence model 
equations integrated. A brief description of this numerical algo- 
rithm and its application to the analysis of basic supersonic 
interactive phenomena was provided in references 9, 10, and 14. 
Details will be provided in this report. 


The application of the above integration options in SCIPVIS in 
analyzing supersonic underexpanded jets exhausting into a super- 
sonic or quiescent external stream will be discussed below. 


4.2 ANALYSIS OF UNDEREXPANDED JET WITH SUPERSONIC EXTERNAL STREAM 


Referring to Figure 9, the analysis of this fully supersonic problem 
(in the absence of embedded subsonic zones) involves utilization 
of the hyperbolic and hyperbolic/parabolic integration options. 

This is accomplished utilizing the techniques summarized below. 
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JET INDUCED BOW SHOCK 


TRANSMITTED 
WAVE 



8M8 M L. 
BOUNDARY 


REFLECTED 

WAVE 


>x,U 


UNDEREXPANDED JET/ SUPERSONIC EXTERNAL STREAM 


FIGURE 9 

Schematic of Mapped and Floating Boundaries Distinguishing 
Various Regions of Supersonic/Subsonic Jet Mixing Solution 


(1) The equations are integrated in mapped rectangular coordi- 
nates utilizing the transformation of equation (26) with 
rL=0 and ru(x) corresponding to the outer edge of the jet 
mixing layer. 


(2) The variation of rjj(x) is obtained by combining the viscous 
and inviscid growth formulations given by equations (33) 
and (35) as follows: 



VIS/ INV 


tan + 



<§ >, 


U 


(36) 
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(3) The position of the inner mixing layer boundary which 

"floats" across the grid points in working its way down 
to the axis is monitored by inspection of the 4> profile 
at each integration step. The hyperbolic (inviscid) 
equations are integrated below this position while the 
hyperbolic/parabolic ( viscous/inviscid) equations are 
integrated above this position (as simply accomplished 
by setting the turbulent viscosity to zero for all grid 
points below the "floating" inner mixing layer boundary*). 


(4) Referring to the insert in Figure 9, properties along the 
outer viscous/inviscid boundary, rjj(x), are obtained via 
solving: 


4 * — 

(a) inner viscous (X vis) and outer inviscid (X inv) 
characteristic compatability relations for the 
pressure, Pu, and flow angle, 0{j; and. 


(b) isentropic streamline relations along the entrained 
streamline, XgL (i.e., H = He, P/pY = const) to 
yield remaining flow properties. 


In original applications of the SCIPVIS model to underexpanded jets 
in a supersonic stream 9 > 10 a single mapping from the axis to the 
fitted jet induced bow shock (Figure 9) was utilized eliminating 
the characteristic jet boundary solution procedure of Step (4) 
above. This was found to be numerically inefficient for a number 
of reasons (to be discussed in the applications section of this 
report). In present applications, the bow shock layer region 
(between ru(x) and the bow shock, tbs( x ) in Figure 9 is analyzed 
utilizing a second mapped domain {here, n = (r - r^(x ) ) / ( tbs( x ) - 
ru(x))} implementating the hyperbolic option to integrate the shock 


* While the hyperbolic/parabolic system of equations could be inte- 
grated across the entire viscous/inviscid jet, the representation 
of turbulence properties in the inviscid core using a high Reynolds 
number, two-equation turbulence model can pose numerical problems 
since the parabolized turbulence equations are not representative 
of the turbulence processes occurring (i.e., the requisite terms 
needed to represent core turbulence, the production of turbulence 
behind shocks and subsequent dissipation, etc. are not included 
in the parabolized turbulence model). Hence, it is most expedi- 
tious to eliminate the solution of the turbulence model equations 
in nonmixing regions rather than to solve an inappropriate system 
of equations. 
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layer equations, or, this region is not analyzed and the external 
supersonic flow is represented via simple one-wave relations such 
as shock-expansion theory (see reference 19). If simple one-wave 
relations are utilized, the outer characteristic compatibility 
relation along X”inV is replaced by a pressure/flow deflection 
rule stating the specific relationship utilized. 


4.3 ANALYSIS OF UNDEREXPANDED JET WITH QUIESCENT EXTERNAL STREAM 


Referring to Figure 10, the analysis of this problem involves utili- 
zation of both the hyperbolic and hyperbolic/parabolic integration 
options required in the fully supersonic problem, plus the parabolic 
or partially-parabolic integration option to analyze the subsonic 
portion of the jet mixing layer. This is accomplished utilizing the 
techniques summarized below. 



UNDEREXPANDED JET/ QUIESCENT EXTERNAL STREAM 


FIGURE 10 
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Schematic of Mapped and Floating Boundaries Distinguishing 
Various Regions of Supersonic/Quiescent Jet Mixing Solution. 



(1) The equations are integrated using the same mapping as 
in the fully supersonic situation described in the 
previous subsection . 

(2) The variation of the outer mixing layer boundary, rjj(x), 
is obtained by combining the viscous growth rule with the 
flow deflection angle at the sonic line position (not the 
outer edge angle as in the supersonic flow case*), viz., 


dr 


U, 


Cr 


dx 


VIS/INV 


= tan0 M=1 + 


U 


/ li ) 

' 3r * 


U 


(37) 


(3) As in the fully supersonic case, the position of the 

"floating" inner mixing layer boundary is used to delin- 
eate between hyperbolic and hyperbolic/parabolic flow 
regions. The hyperbolic/parabolic integration is per- 
formed in the region bounded by the inner mixing layer 
boundary and the jet mixing layer sonic line. 


(4) Above the sonic line, the equations are integrated using 

the parabolic option with the streamwise pressure gradient, 
9P/3x, set equal to zero. The position of the sonic line 
is monitored at each integration step and also "floats" 
through the grid points, ultimately reaching the jet axis 
in the farfield at which point all "steady" wave processes 
terminate . 


(5) Two options are available to represent the radial varia- 
tion of pressure and radial velocity across the subsonic 
portion of the mixing layer which work as follows. 


* For quiescent external flow, the flow angle of the entrained 
streamline can approach -90° and does not represent the 
"effective" inviscid slipstream variation. 
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(a) Parabolic Option 


The pressure variation across the subsonic region 
is neglected. This is a reasonable approximation 
for the quiescent problem at modest pressure ratios; 
it would be a poor approximation for the nearfield 
jet with a subsonic external flow and large boundary 
layer since the streamline curvature induced by the 
"washing away" of the mass defect region of the 
external boundary layer can be quite appreciable 
(see references 3 and 4). In the parabolic option, 
the external pressure is imposed at the subsonic/ 
supersonic matching point (see insert of Figure 10) 
yielding the flow deflection angle via the viscous 
characteristic compatability relation of equation 30. 
Flow properties U, H, and <p at the matching point 
are determined from the parabolic solution. Upward 
integration of continuity with the state constraint 
of equation 34 yields the radial variation of V and 
p and thus yields the outer edge entrainment (radial) 
velocity, Vg, which can be used to upgrade the 
solution in a subsequent iteration as will be 
described below. 


(b) Partially Parabolic Option 

The pressure, P, is assumed at the sonic line yield- 
ing the flow angle, 0, and flow properties as above. 
Upward integration of the normal momentum equation 
yields the radial variation of P. The process is 
iterated upon until the predicted edge pressure is 
equivalent to the imposed external pressure. The 
steps described here are exhibited in Figure 11. 


In performing quiescent external stream calculations, imposition of 
the ambient pressure level along rp and the use of the boundary con- 
dition, Ue = 0 is only a first approximation to the complete inter- 
active problem.* The entrainment velocity., Ve(x) predicted in this 
first approximation accelerates the "nominally" quiescent external 


* The condition Ug = 0 cannot (and should not) be simulated 
numerically since this implies vertical (-90°) angles for 
the entrained streamlines which is physically nonrealistic - 
a nominal value of Ue/Uj ~ 1/100 can be simulated and is 
used to approximate the quiescent state in the first 
approximation . 
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JET MIXING LAYER 



FIGURE 11 

Pressure-Split Analysis of Subsonic Portion 
of Jet Mixing Layer, 


flow. The entrainment induced pressure (Pe(x)) and streamwise 
velocity (Ue(x)) variations (obtained via coupling the jet solution 
with a potential flow solver) are implemented as boundary condi- 
tions along ru for a subsequent jet solution. This jet/potential 
flow iterative coupling approach is analogous to that required 
for analyzing subsonic/transonic external streams. The iterative 
methodology now under development will be discussed in detail in 
Part II of this report since SPLITP is being used to explore 
various coupling procedures. Direct coupling of SCIPVIS with a 
potential flow solver can be accomodated by the present code but 
this capability has not yet been tested. 
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4.4 ANALYSIS OF MACH DISC MIXING REGIONS 


Referring to Figure 12, the mixing process that occurs between 
the inner subsonic flow that has traversed the Mach disc and 
the outer supersonic flow, has strong wake-like characteristics 
(viz. outer to inner pU ratios exceed 10 and can exceed 50 
for highly underexpanded exhausts). The entrainment process 
displaces the outer streamlines inward at a very rapid rate for 
the first few Mach disc radii producing strong interactive 
effects. A simple procedure based on the "porous-sting" approach 
of Dash and Wolf 45 > 6 for the two-phase Mach disc problem has 
been developed for treating this interactive mixing process. 

This simplified approach utilizes a pressure-split "sublayer 
approximation" to eliminate elliptic (saddle-point) behavior 
in the marching process. This approach is applicable for smaller 
radius discs where the extent of the embedded subsonic zone 
behind the disc is small in comparison to the jet width and for 
which the turbulent mixing process dominates the acceleration of 
the inner flow to the supersonic velocities. The results 
obtained are in qualitative accord with the experimental observa- 
tions of Back and Cuffel 47 . The basic features of the Mach disc 
methodology incorporated in SCIPVIS are summarized below. 




FIGURE 12 

Characteristic Features of Mach Disc Mixing Region 
and Subsonic/Supersonic Matching Surface. 
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4.4.1 Disc Location 


A formal triple-point shock fitting calculation is performed at 
each candidate disc location 19 and a disc is dropped when a 
"sting" triple-point criterion is satisfied (i.e., when the 
streamline traversing the triple-point is parallel to the axis; 
locations upstream yield positive angles and locations downstream 
yield negative angles). This yields reasonable disc locations 
for smaller discs but locations somewhat downstream of observed 
locations for discs whose radius exceeds twenty percent of the 
total jet width. If the sting criterion has not been satisfied 
by the time the downrunning shock reaches the jet axis, the shock 
is regularly reflected. In a number of situations, shock reflec- 
tion/disc-formation has been observed 11 * 15 at a position where 
use of the standard triple-point analysis breaks down (i.e., 
where imposition of the normal shock pressure level at the triple- 
point leads to a detached shock solution for the reflected shock) . 
In such situations, the disc is dropped at the observed location 
(which has generally corresponded to the point where the Mach 
number behind the reflected shock at the triple-point first 
becomes sonic) using an "off axis" regular reflection analysis. 
Several of the cases described in reference 15 encountered this 
type of reflection condition and future experiments are planned 
to reveal the details of the specific reflection mechanism occur- 
ring. 


4.4.2 Inner/Outer Matching Boundary 


In cases where a disc is dropped, a grid point, sitting several 
grid points above the viscous sonic line (Figure 12) serves as 
the matching point between the inner subsonic and outer supersonic 
flows. A specified matching Mach number (typically 1.1 to 1.5) 
delineates the grid point at which inner/outer matching is per- 
formed. The grid points above the matching point are integrated 
using the supersonic viscous (hyperbolic/parabolic) or inviscid 
(hyperbolic) option dependent upon the position of the upper Mach 
disc mixing layer boundary (ascertained by inspection of a "dummy" 
inert tracer species profile initially assigned the value of unity 
for flow traversing the disc and zero for all other points). The 
grid points below the matching point are integrated using a 
variant of the partially-parabolic option with the pressure split 
using a sublayer type approximation. The index of the matching 
point changes from step to step in accordance with the monitored 
spatial variation of the sonic line. When the sonic line reaches 
the jet axis, the supersonic equations are employed and the inner/ 
outer matching is no longer required. 
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4.4.3 Inner/Outer Matching 


The following steps are entailed in the inner /outer Mach disc 
matching : 


(1) A value of the pressure gradient, 3P/3x, at the matching 
point 1=1*, is assumed. 


(2) The pressure field is split with the streamwise pressure 
gradient term in the axial momentum equation damped in 
accordance with the following heuristic sublayer approxi- 
mation : 




(Mj* 2 -!) 


for Mj 


> 1 


(38a) 



0 for Mj < 1 


(38b) 


and the pressure, itself, given by: 


P 


I 




Ax 


(38c) 


(i.e., the pressure is uniform across the mixing region 
but the pressure gradient is variable in the transonic 
portion and completely suppressed in the subsonic region). 


(3) The inner grid points are integrated using the partially- 
parabolic algorithm (to be described) yielding the normal 
velocity at the matching surface, V*, and hence the "inner" 
flow deflection angle, 0*. 


(4) The "inner" value of 0* is compared with an "outer" value 
given by a down-running viscous-characteristic compati- 
bility relation (this is analogous to the sonic line 
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matching depicted in Figure 11 with the subsonic and 
supersonic regions inverted). Steps (2) - (4) are 
repeated varying the assumed pressure gradient until 
the inner and outer values of 6* agree to within a 
prescribed tolerance. 


The treatment of larger discs entails accounting for elliptic 
upstream influence effects using iterative methodology in con- 
junction with a downstream saddle point type viscous throat 
constraint 48 analogous to the minimum area throat constraint 
utilized in the Abbett inviscid Mach disc approach 4 * incorporated 
in earlier inviscid versions of the SCIPPY model 2 * 19 . The 
streamwise iterative approach required to treat larger radius 
discs appears to be a straight-forward extension of the present 
small disc methodology. 


4.5 GRID DISTRIBUTION 


4.5.1 Radial Distribution of Grid Points/Embedded Fine Grid 


Referring to Figure 13, the grid points I = 1, 2, ..., IMAX are 
evenly distributed between the jet axis, r = 0, and the outer edge 
of the jet mixing layer, r = rp(x). In the mapped rectangular 
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coordinate system, the transverse coordinate n varies from 0 to 1 
across the jet with the grid spacing, An, given by An = 1/(IMAX-1). 
To provide sufficient grid resolution in the thin nearfield shear 
layer, a grid embedding option is available which subdivides the 
standard grid intervals, An, into a number of subintervals as 
exhibited in Figure 13. 


Referring to Figure 14, the embedded fine grid region is initiated at 
grid point 2 selected to fall one (coarse) grid interval below the 
interval containing the "floating" lower shear layer boundary, 
tl(x). The ratio of fine grid (An) to coarse grid (An) spacing is 
given by : 


= 2 n 


where n is initially user selected. The coarse grid points, 



5 

FIGURE 14 


Details of Embedded Fine Grid. 
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II 


I = I to IMAX, are reindexed in accordance with the relation: 

K = I + (I - I ) • 2 n 


with K varying from I to KMAX * I + (IMAX - I)*2 n as exhibited in 
Figure 14^for n = 2. The coarse grid integration for grid points 
I = 1 to I is performed first yielding properties at £ + AC. The 
radial derivatives at the bounding grid point, I, are evaluated 
using the coarse grid interval, An, and the integration js per- 
formed in the inviscid limit. The interval between I = I and I + 
1 serves as an inviscid buffer zone. 


After the coarse grid integration, the fine grid integration is 
performed from £ to £ + A£ for the grid points K = I + 1 to KMAX 
utilizing a reduced axial step, A£ (chosen to be a fraction of the 
coarse grid step size, A£, in accordance with stability require- 
ments to be discussed below). Note that all properties along K = I 
are known from the coarse grid integration and serve as lower 
boundary conditions for the fine grid integration. At the end of 
the combined fine/coarse grid integration step, the floating lower 
mixing layer edge~ is located and the matching point I is shifted 
if required. As I is moved down, the number of fine grid points 
increases as does the total number of overall grid points. The ratio 
of the shear layer width to the full jet width is maintained to be 
less than l/2 n (i.e., if n=3, the width of the shear layer should 
not exceed fth the overall plume width). If the shear layer width 
ratio exceeds the grid ratio, n is reduced by 1 which decreases the 
number of grid points in the shear layer by a factor of 2 (i.e., 
every other point is discarded). When the shear layer exceeds £ 
the jet width, a single grid interval is used. An additional con- 
straint on n is the total number of grid points dimensioned, viz, 
n must^always be decreased if moving I down leads to the situation 
where I + KMAX > I* where I* is the maximum dimension of I in the 
variable arrays. 


4.5.2 Axial Step-Size Criterion 


The SCIPVIS algorithms are fully-explicit and thus, the allowable 
marching step, A£, is limited by both hyperbolic and parabolic 
stability constraints. The following step-size criterion are 
utilized in the various flow regions of the jet. 
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(a) Hyperbolic Stability Criterion 

In inviscid supersonic flow regions, the step-size is 
limited by the well known CFL condition which requires 
the physical domain of dependence to fall within the 
numerical domain. This criterion is satisfied in a 
non-linear fashion by locating the intersection of the 
A± characteristics with the mapped coordinate grid 
lines, n = const, at each grid point, I, and, deter- 
mining the minimum intersection for all grid points. 

At each grid point I, the maximum allowable step-size 
is given by the relation: 


/ AC ±\ < 

±bAn 

\ I / 

v ' HYP 

tanuij - tan(0±y) I+1 


( 39 ) 


where 0 and y are the flow and Mach angles, and to is 
the mapped grid line angle given by: 


tanto,. = r. 


(I ~ 1) 

(IHAX - 1) 


<v - 


(b) Parabolic Stability Criterion 


In viscous subsonic flow regions, the parabolic "diffusion 
equation" criterion is applied which requires that: 



'5pU0mi n b 2 An 2 


( 40 ) 


(c) Hyperbolic/Parabolic Stability Criterion 

In viscous supersonic flow regions, the hyperbolic and 
parabolic criterion are combined in the following fashion 
at each grid point I : 



HYP/PAR 




( 41 ) 
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This viscous/inviscid step-size was first utilized by Cheng 50 
in the explicit integration of the time-dependent Navier-Stokes 
equations and had proven effective in earlier "viscous- 
characteristic" finite-difference analysis of supersonic mixing 
problems by Dash and coworkers (see, for example, Reference 51). 


4.6 INTERIOR POINT INTEGRATION PROCEDURES 


4.6.1 Generalized Finite-Difference Algorithm 


The following "generalized" two-step algorithm is utilized in all 
regions of the jet flowfield to advance the solution of equation 
(27) from £ to £ + AC at grid point I, for equally space trans- 
verse grid intervals, An: 

Predictor Step: 




(42a) 


(42b) 


where : 


A 


± 



I±1 
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and e is an upwind/alternating one-sided convective difference 
parameter whose operation differs in different flow regions as 
will be described below. The diffusion terms are evaluated 
using a standard central difference operator in all flow regions. 


4.6.2 Supersonic Marching Procedure 


In supersonic viscous or inviscid flow regions, the convective 
difference parameter, e, in the above difference equations is 
varied between zero and one in the predictor and corrector steps 
to provide one-sided, alternating differences for shock- 
capturing. The resultant algorithm is the steady flow equivalent 
of the well known explicit MacCormack algorithm 2 " . The "e" 
sequence is alternated at subsequent marching steps to yield a 
nonpref erential treatment of wave propogation. The pressure 
parameter, a, (see equation 27) is set equal to one in super- 
sonic flow regions, thus including the pressure term in the ep 
conservative array and eliminating it from the source term, G. 

A standard supersonic decode procedure 18 * 19 is used to obtain the 
variables, P, p, U, V, H, <J> , k, e, from the conservation arrays 
at the end of the predictor and corrector steps. 


4.6.3 Upwind Modification of Convective Operator for Scalar 
Variables 


In a number of applications of the supersonic algorithm, oscilla- 
tory behavior of the scalar variables (H, <f>, k and/or e) occurred 
in situations where the mapped grid was significantly skewed with 
respect to the streamline directions and radial gradients of the 
scalar variables were large. In parabolic problems, this unstable 
type of behavior is remedied by the use of an upwind convective 
operator 52 . To provide this remedy in a shock-capturing algorithm, 
the convective difference operator for the scalar components of 
the F vector array can be split into upwind and alternating one- 
sided difference components as follows: 


3F* 

3n 


“ ^ + <> v I § 


alternating upwind 


(43) 


where : 


a 



<J>, k, 


e 




T 
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and F* represents the scalar elements of the F vector array. 

With this splitting, the following convective difference expres- 
sions are utilized in the predictor and corrector steps: 


Predictor Step : 


3F* _ _1_ 
3n ” An 


+ <pv) 


ai ! 

K | (l-e*) a 


(1-e) <pV) I+1 - (l-2e) (pV) x - e(pV> I 




I+1 - (l-2e*) a x - e* 


“i-l | 


Corrector Step : 


3F * = _1_ 
3n = An 


a z |e(pV) I+1 + (l-2e) (pv) x + (e-1) <PV) I _ 1 J 

|e*a I+1 + (l-2e*) a x + (e*-l) | 


+ (0v) x ; e *a I+1 + 


(44a) 


(44b) 


The parameter e is the standard alternating difference parameter, 
whereas e* is the upwind difference parameter defined (in both 
the predictor and corrector steps) by: 


e* = 1 if V/U > a/b 
e* = 0 if V/U < a/b 


(45) 


where a and b are the mapping parameters given by equation 28. 

For the continuity and momentum equations, the standard alternating- 
difference MacCormack algorithm is utilized with no splitting. 

The splitting of the scalar components of the F vector array has 
no effect on the shock-capturing characteristics of the algorithm 
(the scalar variables, a in Equation 43, all change continuously 
across shock-waves). Thus, the nonconservative differencing of these 
variables in the vicinity of the captured shocks poses no problem. 
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This has provided a simple and reliable remedy to the oscillatory 
problems encountered. This approach is recommended as a general 
extension of the conservative MacCormack algorithm in both steady 
and unsteady flow applications, in regions where such oscillatory 
problems are anticipated. In use of the MacCormack algorithm 
with nonconservative flow variables, an analogous practice is the 
upwind solution of the entropy equation in entropy layer regions 
(see references 53 and 54). 


4.6.4 Subsonic Marching Procedure 


For grid points in subsonic regions, the pressure parameter, a, is 
set equal to zero and hence, the streamwise pressure gradient is 
included in the source term array, G, and does not appear in the 
conservation variable arrays E and F. The full system of conser- 
vation equations are integrated from £ to £ + A£ using the two- 
step algorithm of equation 42 in the manner described below. 


(1) Upwind convective differencing (utilizing e* defined by 
equation 45 in 44) is used for all variables in both the 
predictor and corrector steps. 


(2) The pres sure f ield is split with the streamwise pressure 
gradient aP/cix (not 9P/85), imposed and the normal 
pressure gradient, 3P/8rj, evaluated from the known profiles 
at x (predictor) and x+Ax (corrector). 


(3) Dependent variables U, V, H, 4> , k and e are obtained via 
a simple subsonic decode procedure, viz.: 


f = e f /e p ( 46 ) 

where f represents the dependent variables and ep = pU. 

(4) The pressure and density are evaluated from the relation: 


P(£+«), - P(5)j +(!£)«♦§ (HI) AE (47) 

I ' 1 
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( 48 ) 


and the equation of state (listed here for a perfect 
gas) : 

-1 

pU+aOj = - |(U 2 +v 2 ) j ^ 


(5) The conservation variables , ef, are reconstructed in 
accordance with the value of p obtained from equation 
(48) and the decoded variables U, V, H, <p , k, e given 
by equation (46). 


(6) After the predictor and corrector steps have been com- 
pleted, the dependent variables V, P and p are revised 
via a crossflow integration procedure with the variables 
U, H, <J>, k and e held fixed at the corrector values 
established in Step (3). The crossflow integration 
procedure utilizes the continuity and normal momentum 
equations, and the equation of state constraint, and en- 
tails the coupled solution of three nonlinear relations 
at each grid point . Procedures for solving these 
relations are described in detail in Part II of this 
report. The linearized procedure utilized in SCIPVIS 
in situations where the normal pressure variation across 
the subsonic is neglected is described in the next 
subsection. 


(7) The conservation variables, ef, are upgraded in accord- 
ance with the revised values of p and V obtained in 
Step (6). 


4.6.5 Subsonic Continuity Equation Integration 


Let the revised values of p and V at the grid point I,K (point 
C of Figure 15) be given by: 


p = p* + p' 


v = v* + V" 


(49a) 

(49b) 
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FIGURE 15 

Grid Nomenclature for Subsonic Continuity 
Equation Integration. 


where p* and V* designate the values determined via steps (1) to 
(4) of the previous subsection and p' and V' are corrections to 
these values to be determined by the linearized procedure discussed 
below. The mass flux at point C integrating upward from point A is 
given by: 


*C = *A 


+ Ai^ ac + A\p 


AC 


(50) 


where is the flux based upon the unrevised values, viz: 


4 *AC ■ i< p A U A r A J + 'cVcX 1 " 


and 


i(, AC 


is the linearized correction term given by: 




AC 


2 P C U C r C b K An 
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The mass flux at point C integrating streamwise along BC is given 
by: 


*C * *B + A *BC + A *BC 

♦ 

where A\(; BC is based upon the unrevised values, viz: 

4 *BC - - + ‘>C r C J > T C* 5 

+ 5^ p B U B r B J + P C U C r C ^ [ax] 45 

L J BC 

with the grid line inclination, dr/dx given by: 

[S] bc ' 5 [ r L + "< r U- r L>] K + I [ r L + "< r D- r L>] K . x 


and the 


linearized correction term, 


A^ bc , given by: 



Equating the values of given by equations (50) and (51) 
the linear p", V" relation listed below: 


A 1 P C 


A 2 V C - A ^ 


BC A ^AC 


+ * B ~ ♦ *) 


( 51 ) 


yields 

(52) 


49 



where : 


A 1 * r C U C 


b K 4n + ( V d) c ' [«*] 


BC 


and : 


A 


2 


p B r B J + 


* J 

p c r c 


Replacing Vc in equation (52) with the linearized equation of 
state constraint given by: 


v c - 


H c -i(U c 2 +V* 2 ) 

P& 


Pr- - 


yields : 

- _ 2 a 4c~ a ^ac + ^b-^a 

P c (A 1+ A 2 /A 3 )A£ 


(53) 


(54) 


The subsonic continuity integration sweeps upward from the axis, 
for the Mach disc subsonic region, and from the jet sonic line, 
for the jet mixing layer subsonic region. Both integrations are 
initialized with p' = V' = 0. 


It should be noted that the streamwise continuity relation given 
by equation (51) does not involve the V distribution at the pre- 
vious streamwise step as is appropriate to parabolic situations 55 ’ 56 . 
This implies that the V distribution is evaluated at the midpoint 
of the integration step (i.e., £ + A£/2), an assumption generally 
invoked in present generation parabolic algorithms. The use of 
the normal momentum equation to yield the provisional values of V 
in situations where the normal pressure variation is neglected 
serves as a "redundant" stabilizing factor in the overall algo- 
rithm and essentially implies that the normal velocity component 
is convected along the flow streamlines until corrected by the 
linearized continuity relations of this subsection. The basic 
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assumption is that the continuity induced V" correction to the 
value predicted by the normal momentum equation is small, which 
is generally the case except in the initial jet region, and in 
regions at the end of shock cells where the flow angle changes 
discontinuously in the vicinity of the jet sonic line. In these 
regions, several iterations may be required accompanied by an 
underrelaxation process to stabilize the V" change from iteration 
to iteration. 


4.7 BOUNDARY POINT PROCEDURES 


The following boundary point procedures are incorporated in the 
SCIPVIS model: 


(1) shock point 

(2) wall point 

(3) axis calculation 

(4) jet mixing layer outer edge 

(a) supersonic external stream 

(b) subsonic external stream 

(5) subsonic/supersonic matching point 

(a) jet mixing layer 

(b) Mach disc mixing layer 


The shock point and wall point calculations are performed in the 
inviscid limit* utilizing predictor/corrector characteristic based 
procedures which are described in detail in references 57 - 59. 

The remaining boundary point procedures will be described below. 


* The shock point calculation is used to analyze the jet induced 
bow shock depicted in Figure 9; the wall point calculation is 
utilized in analyzing ducted mixing problems (to be described 
in Section 5) with the wall region treated in an inviscid fashion 
(i.e., the mixing region does not extend to the duct walls and 
wall boundary layer effects are neglected). 
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4.7.1 Axis Calculation 

At the jet axis (or plane of symmetry), the conservation equations 
take the following limiting form: 


3E 

3X 


(1+J)pf || 


+ G. = (1+J) 


y 


eff 


3 2 f 

3r 2 


(55) 


where : 

3V = b ( V 2~ V l) 
3r Aq 


and : 



where subscripts 1 and 2 denote the grid points at the_axis and^ 
one grid interval above the axis. The vector arrays, E, f and Gf 
are defined in Section 3.5. A two-step predictor/corrector pro- 
cedure is used to integrate equation (55) in both supersonic and 
subsonic regions (behind Mach discs) with the pressure parameter 
a set equal to 1 or 0 and the axial pressure gradient predicted 
or prescribed accordingly. 


The supersonic axis solution procedure is modified beyond the 
position of the first reflected shock (or Mach disc sonic line) to 
avoid oscillatory behavior at subsequent shock reflection points. 
The modification effectively reflects waves off of a sting whose 
radius corresponds to one grid interval. This eliminates the 
highly nonlinear axisymmetric wave strengthening that occurs in 
the vicinity of the axis that cannot be resolved using standard 
grid spacing. The modified axis procedure works as follows: 
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(1) The axial pressure gradient at grid point 2 is deter- 
mined utilizing a X” viscous-characteristic relation 
(see equation 30) with the flow deflection angle at 
point 2 set equal to zero (i.e., the flow is enforced 
to be parallel to the axis). 


(2) The conservation equations at grid points 1 and 2 are 
solved with a = 0 imposing the axial pressure gradient 
determined from the X“ relation of step (1) and 
neglecting the normal pressure variation between 1 and 
2. 


(3) The conservation variables at grid points 1 and 2 are 
decoded using the simple (pressure prescribed) decode 
relations of equation 46; V is set equal to zero at 
both grid points and the density is determined from 
the state relation (equation 34). 


4.7.2 Jet Mixing Layer Outer Edge Calculation 


For a supersonic external stream (Figure 9), properties along the 
outer edge boundary, ru(x), are determined using characteristic 
relations as exhibited in the insert of Figure 9 and summarized 
in Section 4.2. The calculation is analogous to a standard con- 
tact surface calculation (see references 19, 57 - 59) modified 
to account for jet entrainment across the boundary. The following 
sequence of operations is entailed in performing the supersonic 
outer edge calculation in marching from £ to £ + A£: 


(1) The outer edge is locating utilizing equation 36, which 
takes the following difference form for I = IMAX: 


Predictor : 


*U = r U (5) + 


tane^S) + CryCOMD f (Q Ar) 


>A5 


Corrector : 


ru (5+A?) - f 


r u (5)+? u + 


tanSy + Cr y b 


( f r f i-i) 

f T An 


A£ 


(56a) 


(56b) 


53 


( 2 ) 


The "feet" of the X VIS> and Xsl characteristics 
(i.e., the intersection of these characteristics with 
the initial data line at £) are determined via com- 
bining the characteristic geometry relations given 
below with linear property variation relations along 
the initial data line. The characteristic geometry 
relations are given by: 


Predictor : 


r*(£) = r D + A*U)A£ (57a) 


Corrector : 


r*(£) = r v a+ AC) + | X*(0+X(S+AO A? (57b) 


and are solved concurrently with interpolation relations 
for X* in both the predictor and corrector steps to yield 
a noniterative location procedure (see references 19 and 
59 for details). In the above relations, X* represents 
A + , X - or Xsl where X designates the slope of the charact- 
eristic line (viz., A± = tan(0±y); Xsl = tan0). Properties 
at the foot locations, r*(£), are obtained via linear 
interpolation. The viscous source terms in the X+vis com- 
patability relation are evaluated via standard central 
difference relations at the grid points and interpolated 
at the foot location. 


(3) The pressure, P, and flow deflection angle, 0, are eval- 
uated via the A + and X" compatibility relations listed 
below. 

Predictor : 

A + (5) fcn (P/P + ) + 0-0 + = B + (£)A5 
A"(C) tn (P/P") - 0+0“ = B“(£)A£ 


(58a) 

(58b) 
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Corrector : 


j |£*(0+a + (?+A 0J *n(P/P) + 

y £a"(C)+A'(C+AC)] *n(P/P + ) - 
In the above relations, 


•4. 

o-fr 


e+e 


y [b 4- (C)+B + (C+AC) J AC 
y [b _ (C)+b'(C+AC)] AC 


(59a) 


(59b) 


„± _ siny cosy 

a - Y 


and: 


+ 


B” 


-Jsin6 

r 



siny 

cos( 0+y) 


+ + + 

In the predictor step, A and B , and, the properties 0 , 

P- in equation (58) are evaluated at the "foot" points r- 
given by equation (57). In the corrector step. A* and B* 
are averaged across the characteristic, and, the ~ nota- 
tion for the values at the C station indicates that these 
values have been updated in accordance with the shift in 
the characeristic foot position between the predictor and 
corrector steps. „ The same nomenclature applies to the up- 
dated values of 0* and P ± at the C station. 


(4) The total enthalpy and density are evaluated along the 
streamline Ysl v * a t ^ ie relations: 


dH 

dX 


0 


and : 


d iZlsll = 0 

dX 
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(5) The species parameter, <J>, the turbulent kinetic energy, 
k, and the turbulent parameters e or W are set equal 

to zero along the outer boundary. 

(6) The velocity, Q=(u 2 +V 2 ) is determined from the state 
relation (equation 34) and the U and V components are 
evaluated from the relations: 


U = Qcos0 
V = Qsin8 


If the external flow is simply represented via the use of 
pressure/flow-deflection relations, equations 58b and 59b are 
replaced by these relations, and the entropy level is assumed 
constant in the external flow. 


For subsonic or quiescent external flow (Figure 10), the axial 
pressure variation is prescribed along the outer jet boundary, 
rjj(x), along with the axial velocity variation, Uu( x ) . The 
total, enthalpy is taken to be constant and <t>, k and e (or W) are 
all set to zero. The radial (entrainment) velocity variation, 
Ve(x), is determined via integration of the continuity equation 
aposteriori, as described in Section 4.6.5. The coupling of the 
subsonic and supersonic portions of the jet at the viscous sonic 
line (Figure 10) entails the use of characteristic methodology 
as described below. 


4.7.3 Subsonic/Supersonic Matching Point 


Referring to the insert of Figure 10, subsonic/supersonic coup- 
ling at the matching point, I*, entails stipulation of the 
pressure at the matching point and determination of the flow 
angle, 0*, via a A + viscous characteristic relation. For the 
quiescent external flow problem in a mildly underexpanded situa- 
tion, the external pressure is imposed at the matching point. 

Use of the subsonic marching procedure described in Section 4.6.4 
yields the dependent variables U, V, IF, <p , k and e at the match- 
ing point. The normal velocity, V, is then revised in accordance 
with the relation: 

V* = U tan0* 
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and the density, p*, is determined from the state relation. The 
subsonic continuity integration procedure (Section 4.6.5) is 
initiated at the I* point and works upward to the jet outer edge. 


For subsonic external flow, the pressure at the matching point is 
assumed but the axial pressure gradient at the outer edge (SP/dx)jj 
is imposed in performing the subsonic marching procedure along 
I* (viz., the pressure field is split in the standard fashion). 

In the first iteration, the pressure at the matching point, P*, 
is taken to be given by: 


P*(5+A£) = P*(£) 



A? 


The iteration proceeds until the pressure at the jet outer edge 
determined from the coupled solution of the continuity and normal 
momentum equations agrees with the prescribed outer edge pressure 
This iterative methodology has been schematized in Figure 11. 


An analogous procedure is performed at the subsonic/supersonic 
matching point in the Mach disc mixing layer (Figure 12). The 
pressure splitting differs here since the upstream influence 
effects can produce departure solutions unless the axial pressure 
gradient is supressed using a sublayer type approximation. In 
the jet mixing layer, the upstream influence effect is already 
incorporated via stipulation of the axial pressure gradient in 
accordance with an external potential flow solution arrived at 
apriori. The details of the Mach disc matching point solution 
have been provided in Section 4.4.2. 


4.8 INITIALIZATION PROCEDURE 


The initialization procedure utilized is directly analogous to 
that employed in the previous overlaid viscous/inviscid jet 
models 1 * . The initialization entails: 


(1) determining the initial jet inclination angle in the 

inviscid limit to achieve balanced pressure across the 
jet slipstream (see references 18 and 59 for details); 
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(2) determining initial mean flow shear layer profiles 
centered about the jet slipstream location a small 
distance downstream of the nozzle exit plane (see 
references 1 and 2 for details); and, 


(3) determining initial turbulence parameter (k and e 
or k and W) profiles based on an equilibrium tur- 
bulence hypothesis and the use of a mixing length 
turbulence model approximation (see references 1, 2 
and 30 for details) . 
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5.0 APPLICATIONS OF SCIPVIS MODEL 


5.1 WAVE/SHEAR LAYER INTERACTIONS 


To illustrate some of the fundamental interactive phenomena that 
can occur in the nearfield jet shear layer, a series of calcula- 
tions were performed for a supersonic shear layer in a duct. The 
wall region of the duct was treated inviscidly and the duct wall 
was turned to generate shock and expansion waves. All calculations 
were performed using the kel turbulence model and employed thirty 
grid points equally spaced between the duct walls. 


SUPERSONIC/SUPERSONIC SHEAR LAYER 



FIGURE 16 

Predicted Wave Pattern Induced by Supersonic Shear Layer 
in a Parallel Wall Duct. 


5.1.1 Waves Generated by Supersonic Shear Layer 


In the first calculation, the duct walls remained parallel. Ex- 
hibited in Figure 16 are the perturbation pressure waves produced 
by the mixing of uniform Mach 2.4 and Mach 4.8 streams, both 
initially at T = 1000 °K and P = 1 atm. Figure 17 depicts the 
mixing induced pressure variations along the duct walls. As dis- 
cussed earlier, in high speed mixing problems such as this, 
significant viscous (turbulent) dissipation occurs. This produces 
a peak temperature in the mixing layer and hence, a hot central 
region of lower density. Streamlines are diverted around this hot 
central region producing compression waves, As the mixing layer 
grows, the extent of the hot central region increases and the 
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FIGURE 17 

Mixing Induced Pressure Variations Along Upper 
and Lower Duct Walls. 


pressure must increase correspondingly. This hot central zone 
has a local mass defect producing a positive displacement effect 
akin to that of a boundary layer. The displacement thickness 
increases with the growth of the mixing region. In supersonic 
combustion problems^ 1 , or in situations with larger velocity 
ratios or higher Mach numbers, the heat produced will be much 
greater and shock waves will be generated. 


The initial wave fronts are clearly displayed in the temperature 
profile at x = 1 (Figure 18) which shows a sharp distinction 
between the mixing region and the perturbed outer flow. The wave 
fronts in the initial region (0 < x < 2) are also readily dis- 
played by the normal velocity profiles depicted in Figure 19. 
Inspection of the profiles indicates that the normal velocity is 
essentially zero in the center of the mixing layer, positive 
above this point and negative below it. Hence, the streamlines 
are seen to be diverted about this central point where the temp- 
erature peaks. Note that maximum and minimum values of the normal 
velocity occur near the edges of the mixing layer. Outside the 
mixing layer, the streamlines are deflected due to the mixing 
induced compression systems. Consider the profile at x = .5. 

The upper and lower wave fronts occur at Y = .75 and .3. At x *= 1, 
the upper wave front has reached the wall while the lower wave 
front is at y = .2. At x = 2, the lower wave front has just 
reached the wall. 
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FIGURE 20 

Temperature Profiles at Several Axial Stations. 


Temperature profiles are displayed in Figure 20 at various axial 
stations. These profiles display the peak value in the center 
of the mixing layer and clearly distinguish the mixing layer and 
wave zones. At x = 2, the temperature reaches the appropriate 
peak value of about 1300 °K but the profile is sharp at the peak 
point and not yet fully developed. Here, only 5 grid points are 
within the mixing zone and thus the resolution is quite poor. 

The lower wave front is seen to reach the wall at this station. 

At x = 4, the profile starts to take on a rounded fully-developed 
appearance. The peak values at later stations go beyond 1300 °K 
due to the continual compression waves in the flow. Profiles of 
the mass fraction parameter, 4> , at various axial stations are 
exhibited in Figure 21 and are essentially identical to the axial 
velocity profiles depicted in Figure 22 (the perturbation pressure 
waves have a negligible influence on the axial velocity and the 
prediction was performed with a turbulent Prandtl number of unity). 


Profiles of the turbulent kinetic energy and dissipation rate are 
shown in Figures 23 and 24 respectively. The initial level of k 
as given in a mixing length initialization is seen to be under- 
estimated. Since the width of the initial mixing layer is quite 
small, this should have a negligible influence on the solution. 
The peak value of turbulent kinetic energy stabilizes at x > 4. 
The profile at x = 2 overshoots this value indicating that the 
solution here is not yet fully developed. The peak value rep- 
resents maximum velocity fluctuations of about 12.7%. 
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The predictions indicate that the mean flow profiles do not 
become fully-developed until x > 2 while the turbulence profiles 
do not until x > 4. This is due to lack of grid resolution in 
the shear layer. Use of the grid embedding procedure described 
in Section.4 could have been employed to substantially increase the 
resolution , 


5.1.2 Expansion Fan/Shear Layer Interaction 


This calculation has the same initial conditions and initial duct 
geometry as the above case. At x = 2, the upper wall is turned 
up by .2 radians generating a downrunning expansion fan which in- 
teracts with the shear layer. The interaction process involves 
the partial transmission of the expansion fan through the shear 
layer; the partial reflection of the fan from the shear layer; and, 
the alteration of mean flow and turbulence variables in the shear 
layer due to the expansion waves. 


EXPANSION FAN/SHEAR LAYER INTERACTION 



FIGURE 25 

Predicted Flow Pattern for Interaction of Expansion Fan 
with a Supersonic Shear Layer in a Duct . 


The predicted flow pattern is depicted in Figure 25 and wall pres- 
sure variations are shown in Figure 26, The pressure along the 
upper wall is perturbed at x = 1 due to the mixing induced weak 
compressions and drops abruptly at the expansion corner (x = 2) by 
a factor of about 2.3. The pressure then remains constant along 
the upper wall until the wave front of the reflected fan reaches 
the wall at x ~ 3,5. At x * 8, the interaction of the reflected 
fan with the upper wall appears completed. The pressure along the 
lower wall is perturbed at x ~ 2 by the mixing induced compression 
with the front of the transmitted expansion wave reaching the lower 
wall at x ~ 5. 
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FIGURE 26 

Pressure Variations Along Upper and Lower 
Duct Walls. 


The sequence of events is readily portrayed by inspection of the 
temperature profiles (Figure 27) at various axial stations, At 
x = 2, just downstream of the corner expansion, the wall temp- 
erature decreases to 800 °K behind the expansion fan. At x = 

2.5, the expansion fan is between .7 < y < .95 and is about to 
intersect the shear layer. The reflected wave front is at y - 
.95 but the strength of the initial reflected waves are quite 
small. At x = 3 , the front of the expansion fan has traversed 
half the shear layer and has reduced the peak temperature by 
about 50 °K. At x = 3.5, the front of the expansion fan has tra- 
versed the entire shear layer and the peak temperature is now 
1175 °K (reduced by 125 °K). At x = 4.0, the expansion fan is 
between .2 < y < .65. Note that the peak temperature has dropped 
to 1060 °K (by 240 °K). The reflected waves generated in the 
interaction process are now starting to become appreciable as 
evidenced by temperature levels below 800 °K at y ~ .89. At x = 

6, the expansion fan front has reflected off the lower wall and 
is located at y ~ .1. The peak temperature level has dropped to 
975 ok. Temperatures at the upper wall are now lowered to 750 °K 
resulting from the wave reflected from the shear layer. At x = 8 , 
the expansion fan front is at y ~ .5., while at x = 10, it is at 
y ~ .8. 
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5.1.3 Shock/Shear Layer Interaction 


This calculation is analogous to that above with the upper wall 
turned down at x = 2 by .2 radians. This generates a downrunning 
shock wave which interacts with the shear layer. The interaction 
process involves the partial transmission of the shock through the 
shear layer; the partial reflection of the shock from the shear 
layer; and the alteration of mean flow and turbulence variables in 
the shear layer due to the shock wave. 


SHOCK /SHEAR LAYER INTERACTION 


Y 



FIGURE 28 

Predicted Flow Pattern for Interaction of Shock Wave 
with a Supersonic Shear Layer in a Duct , 


The predicted flow pattern is depicted in Figure 28 and wall pres- 
sure variations are shown in Figure 29. The pressure perturbations 
induced by the initial mixing region are not distinguisable on this 
scale. The pressure along the upper wall increases by a factor of 
about 2 at the compression corner (x = 2), It then remains constant 
until reflected compressions start impinging at x ~ 3 . The pressure 
along the lower wall remains constant until the downrunning shock 
generated at the compression corner impinges at x ~ 4.2. The re- 
flected compression waves generated by the interaction of the shock 
with the shear layer reflect from the upper wall as compressions and 
coalesce to form a shock in propagating to the lower wall. This 
second shock reaches the lower wall at x - 4.6. The two shocks 
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FIGURE 29 

Pressure Variations Along Upper and Lower Duct 
Walls. 


reflect off the lower wall and intersect in propagating back to the 
upper wall. They reach the upper wall at x ~ 5.2 at which point 
the calculation is terminated since the flow behind this reflection 
point become subsonic. 


This sequence of events is readily monitored by inspection of the 
pressure profiles at various axial stations (Figure 30). At x = 2 , 
gust downstream of the compression corner, the wall pressure is 
seen to have increased by a factor of 2. At x = 2.5, the captured 
shock front has traversed about one-fourth the shear layer. The 
interaction process is clearly displayed in this figure. The pres- 
sure level is seen to be constant (P ~ 2) downstream of the shock 
outside of the shear layer but increases in the shear layer due to 
the reflected compression waves generated. At x = 3 , the shock has 
traversed the entire shear layer. The positions of the captured 
shock and the reflected compression waves are indicated. At x = 
3.5, the front of the captured shock is at y - .2 and the compres- 
sions reflected from the upper wall are between .4 < y < .6. 
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At x = 4, the primary shock is about to impinge on the lower wall 
while the compressions have focused to yield a secondary embedded 
shock located between .17 < y < .22. At x = 4.5, the primary 
shock has reflected bringing the wall pressure to 5.25. The junc- 
ture of the downrunning secondary shock and uprunning reflected 
primary shock is at y ~ .025 (where the pressure peaks at 6). At 
x = 5, the reflected primary shock is between . 15 < y < .24 and 
is somewhat diffuse. The reflected secondary shock is between .05 
< y < .1. 


These same events are portrayed in the temperature profiles given 
in Figure 31. These profiles can be more informative then the 
pressure profiles since the shear layer region can be readily iden- 
tified which helps clarify interactive events. The shear layer 
boundaries are labeled by 4<i and in all temperature profiles. 
Note that after the primary shock traverses the shear layer at x = 
3, the peak temperature increases to 1700 °K. After the secondary 
shock traverses at x = 4 , the temperature increases to about 1750. 
With the passage of the primary shock reflected from the lower wall 
at x = 5, the temperature increases to 2150. 


Further details of these fundamental interactive calculations are 
provided in reference 9. 


5.2 LOW SPEED PARABOLIC JET MIXING 


The SCIPVIS model was developed as an extension of the inviscid 
SCIPPY model 18 ’ 19 . Hence, while its performance in analyzing hyper- 
bolic wave processes is quite well established, its performance in 
analyzing parabolic turbulent jet mixing processes requires detail- 
ed assessment. This assessment was performed by comparing SCIPVIS 
"parabolic mode" predictions with those of the well established 
parabolic jet mixing model, BOAT 1 (which has been used in a number 
of turbulence model evaluation studies 30 ’ 90 ’ 61 encompassing a wide 
body of balanced pressure imcompressible and compressible free jet 
mixing problems) and the new jet mixing model, SPLITP (whose 
features will be described in Part II). 


5.2.1 Low Speed Jet into Moving Stream, Uj/Ug = 2 

Comparisons of SCIPVIS predictions with those of BOAT and SPLITP 
for a simple incompressible jet with a velocity ratio Uj/Ug = 2 
are exhibited in Figures 32 and 33, The kel model was utilized 
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FIGURE 32 

Comparison of BOAT, SPLITP and SCIPVIS Parabolic Predictions 
for Centerline Decay in Low Speed Jet Mixing. 



Comparison of BOAT, SPLITP and SCIPVIS Parabolic Predictions 
for Turbulent Kinetic Energy Variation Along Centerline in 
Low Speed Jet Mixing. 
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for all predictions and 41 grid points were employed. The predic- 
tions for the centerline variation of the species mass fraction 
parameter, shown in Figure 32, (initially, <j> = 1 in the jet and 0 
in the external stream) and the turbulent kinetic energy, shown in 
Figure 33, are all in close agreement. Differences are attributed 
to the different grid distributions utilized (i.e., BOAT employs a 
mapped stream-function distribution which concentrates all the grid 
points in the thin nearfield shear layer region while SCIPVIS 
(without the grid embedding option) evenly distributes grid points 
from the jet axis to the outer mixing layer edge and thus has 
coarse shear layer resolution). This does not, however, signifi- 
cantly effect the quality of the nearfield prediction as exhibited 
in Figure 34. Here, SCIPVIS with only 8 of its 41 grid points in 
the shear layer (at the axial station x/rj = 6) predicts profile 
shapes and peak turbulence levels which are. in accord with those 
predicted by BOAT with all 41 points distributed across the shear 
layer. 



Comparison of Nearfield Species Parameter and Turbulent 
Kinetic Energy Profile Predictions. 
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5.2.2 Low Speed Jet into Still Air 


In low speed jet mixing problems, the performance of the kW 35,1,2 
and ke 31 turbulence models have sometimes proven comparable 30 , as 
exhibited in Figure 35 for an in compressible jet into still air. 



FIGURE 35 


Comparison Between Predictions and Data for an Incompressible 
Jet with Specified Initial Velocity Profile. 


Here, the predictions were initialized several jet radii downstream 
(i.e., at x/rj ~ 2) using the mean flow measurements of Maestrello 
and McDaid 62 and a mixing length approach to initialize the turbu- 
lent variables. The agreement with the centerline velocity decay 
data is quite good and the two predictions are virtually indistin- 
guishable from each other. However, the kW model has been found 
to poorly predict nearfield mixing rates (see Figure 6) and thus 
cannot correctly predict jet core lengths unless initial mean flow 
data is prescribed some radii downstream of the nozzle (viz. , 
operation of the kW model with a top hat profile at the nozzle exit 
plane predicts overly long jet core lengths). However, the ke 
models accurately calculate incompressible 2-D shear layers and, 
because of the use of an axisymmetric correction, also calculate 
the developed region of jets with good accuracy (see reference 30). 
Hence, these models are recommended for SCIPVIS predictions of low 
speed jet mixing. 
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5.3 SUPERSONIC BALANCED PRESSURE JET MIXING 


SCIPVIS can perform supersonic balanced pressure jet mixing prob- 
lems in either a parabolic or fully-coupled viscous/inviscid mode. 
Unless the weak waves induced by the mixing process are of interest, 
the results using either of these modes are virtually identical. 


U/Uj 



FIGURE 36 

Comparison of Fully-Coupled SCIPVIS and Parabolic Boat 
BOAT Predictions for Balanced Pressure Mach 2.2 Jet 
into Still Air. 


Figure 36 compares a fully-coupled SCIPVIS prediction for a balanced 
pressure Mach 2.2 jet into still air with a parabolic BOAT predic- 
tion. The ke2 model was implemented for these runs. A prediction 
made running SCIPVIS in the parabolic mode was virtually identical 
to the fully-coupled prediction except for the wavelets in the core 
region and thus, is not exhibited. The wavelets are physical (see 
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the measurements of McLaughlin, et al. 63 which exhibits the meas- 
ured centerline Mach number variation for a Mach 2 jet into still 
air) and are generated by the interaction mechanism discussed 
previously . 


5.3.1 Analysis of Eggers Data for Mach 2.2 Jet into Still Air 


The performance of the three turbulence models in SCIPVIS (para- 
bolic mode) in predicting the decay in centerline velocity for the 
Mach 2.2 jet is exhibited in Figure 37 and compared with the data 
of Eggers 64 . Velocity profiles at x/rj = 8 and 27 are compared 
with data in Figure 38. The comparisons here indicate that the 
kW turbulence model does the best overall job, the ke2 model 
mixes too fast, and the ke2,cc model mixes somewhat too slowly. 



FIGURE 37 

Comparison Between Predicted and Measured Centerline Velocity 
Decay for a Mach 2.2 Jet into Still Air, 
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• Eggers Data 




FIGURE 38 

Comparison Between Predicted and Measured Radial Velocity 
Profiles for a Mach 2.2 Jet into Still Air, 


The available data for balanced pressure, supersonic jet mixing is 
quite sparse and does not permit adequately assessing the general- 
ized performance of the above turbulence models. An experimental 
program has recently been completed at CALSPAN^ 5 which signifi- 
cantly extends this data base (viz., data has been obtained for 
Mach 2, 3, and 4 jets into quiescent and supersonic external 
streams with density ratios varied by varying both temperature and 
composition). A program (supported by the Air Force Rocket Propul- 
sion Laboratory) is now in progress by Pergament and coworkers at 
Science Applications, Inc. to use this new data in assessing the 
above turbulence models. 


5.3.2 Analysis of Beach Data for Hg Jet, Free and Ducted Mixing 

Calculations were performed to simulate experiments of Beach 66 at 
NASA/LRC for a cold Mach 2 H 2 jet exhausting into a hot Mach 1.9 
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N 2 stream. Test cases 2 and 5 of reference 68 were analyzed for 
nonreacting, balanced pressure mixing. The injector used had a 
substantial lip thickness (~ 4 of the jet radius) and the jet ex- 
haust was conical (half-angle of 5 degrees). A large external 
boundary layer (6 ~ rj) was present. In test case 2, the mixing 
was free with the N 2 stream extending out to 6.87 radii. In test 
case 5, the mixing was ducted with the duct radius equalling 6.87 
jet radii. 


PREDICTED FLOW STRUCTURE FOR 



Predicted Flow Pattern for Mach 2 H 2 Jet Exhausting into 
M = 1.9 N 2 Stream. 


The flow pattern predicted by SCIPVIS is shown in Figure 39. A 
multi-cell wave structure occurs in the jet primarily induced by 
the conical exhaust conditions. The waves in the jet are partially 
transmitted through the mixing layer and are reflected off the 
interface (treated inviscidly) separating the N 2 stream from the 
quiescent environment. These reflected waves reach the axis at 
x/rj ~ 20 and the second cycle of these waves reaches the axis at 
x/rj ~ 40. This is exhibited in Figure 39 and is also made evident 
by the predicted centerline pressure variation shown in Figure 40. 
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FIGURE 40 

Centerline Pressure Variation. 


The initial pitot pressures for cases 2 and 5 are identical and 
are shown in Figure 41. The darkened and open data points of 
Beach correspond to measurements on either side of the jet center- 
line. This data was fitted by the dashed curve to initialize the 
calculation. The initial total enthalpy and species variations 
were taken to be step profiles with the step occurring at the 
midpoint of the injector lip. With these profiles stipulated and 
the pressure assumed to be ambient, the profile for Mach number 
shown in Figure 41 was obtained. The initial turbulent kinetic 
energy profile was obtained using the dual length scale mixing 
length formulation described in Reference 1. The Mach number dis- 
tribution along the lip face (between .65 < r/rj < .90) is subsonic 
with the minimum Mach number of .2. The fully-coupled option in 
SCIPVIS does not yet treat off-axis embedded subsonic zones such as 
this; hence, the initial pitot pressure profile was modified as 
shown in Figure 41 to boost the flow "artificially" to supersonic 
velocities (i.e., a constant pressure level of 2.15 atm. was used 
along the lip face). This boosted the Mach numbers to 1.1 across 
the face. 
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Modified and Unmodified Initial Profiles for H 2 /N 2 Jet 
Mixing Case. 



FIGURE 42 

Comparison of Predicted Pitot-Pressure Variation Along 
Jet Centerline with Beach Data; Predictions Performed 
with Actual and Modified Initial Profiles using Para- 
bolic and Fully-Coupled Run Options. 
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To ascertain the sensitivities to altering the initial profiles 
and to assess the differences between the fully-coupled and 
parabolic (constant pressure) solutions, the following three 
calculations were performed : 


A. Parabolic SCIPVIS run - unmodified initial profiles 

B. Parabolic SCIPVIS run - modified initial profiles 

C. Fully-Coupled SCIPVIS run - modified initial profiles 


The predicted variations of pitot-pressure along the jet center- 
line are compared with the Beach data in Figure 42. The parabolic 
run with the unmodified initial profiles agrees quite well with 
the data. The parabolic and fully-coupled runs with the modified 
profiles initially mix more slowly since the initial turbulence 
levels are substantially reduced and the low velocities in the lip 
region (which enhance the mixing) have been artifically boosted to 
a Mach 1.1 level. The fully-coupled and parabolic predictions ex- 
hibit the same rate of mixing. We conclude that the fully-coupled 
prediction would agree with the centerline decay data using the 
correct initial conditions, To perform this calculation, we 
require a means of dealing with the embedded subsonic region. 
Pressure-splitting type subsonic marching can be gainfully employed 
for this purpose. 

In Case 5 of reference 66, the mixing is ducted. SCIPVIS was run 
in the fully-coupled mode with the modified initial conditions. 

To account for wall boundary layer effects (SCIPVIS implements 
inviscid wall conditions) the duct wall (r/rj = 6.875) was modi- 
fied by a linear 6*(x) distribution (6*(x) = .003875 x/rj) so that 
the duct wall varied from 6.875 to 6.5 in 100 jet radii. The cal- 
culated pitot-pressures at the end of the duct (x/rj = 100) are 
compared with data and with GENMIX parabolic model predictions in 
Figure 43. The GENMIX and SCIPVIS predictions are virtually iden- 
tical except in the vicinity of the wall where SCIPVIS does not 
calculate the boundary layer. Both agree quite well with the data. 
(The GENMIX predictions were extracted from reference 66 and 
utilized the same ke turbulence model formulation as the SCIPVIS 
prediction). Further details of both these calculations are avail- 
able in reference 9. 
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FIGURE 43 

Comparison of Pt Profile at x/rj = 96 with Beach Data 
and GENMIX Predictions. 


5.4 UNDEREXPANDED JETS INTO SUPERSONIC EXTERNAL STREAM 


5.4.1 Mach 3 Jet into Mach 3 External Stream; Axis to Bow Shock 
Mapping 


A series of calculations were performed for underexpanded super- 
sonic jets exhausting into supersonic external streams to assess 
the SCIPVIS code's ability to analyze plumes with strong discontin- 
uities. The first calculation was for a Mach 3 jet exhausting 
into a Mach 3 external stream at a static pressure ratio of 10/1. 
The jet temperature was 1800 °K and the external stream temperature 
was 300 °K. The calculation was performed using a single mapped 
domain extending from the axis to a fitted bow shock - 50 radial 
grid points were implemented. The calculation was initiated using 
an overlaid composite solution at x/rj = .5. The predicted flow 
structure for the first shock cell is shown in Figure 44. The 
positions of the embedded shear layer boundaries, i|)j and \p 2 , the 
viscous dividing streamline, 4> 0 , and the embedded barrel and re- 
flected shocks are all indicated. The orientation of the mapped 
coordinate lines n=0, .1, ... .9 and 1.0 is also indicated. 

These lines intersect the viscous dividing streamline at moderate 
angles (i.e., at x/rj - 14, the intersection angle is about 15°) 
and thus the mapping is not ideal. An ideal mapping would closely 
align the coordinate lines with the streamlines to maximize the 
local Courant numbers and to minimize convective errors). 


83 




6 

K * 10* 


FIGURE 45 

Mass Fraction and Turbulent Kinetic Energy Profiles at x/rj = 
12.4 Exhibiting Extent of Mixing Layer. 
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Profiles of the mass fraction parameter, <j>, and the turbulent 
kinetic energy, k, at x/rj = 12.4 are depicted in Figure 45 and 
clearly delineate the boundaries of the mixing region. Less than 
10 grid points fall within the mixing region at this station yet 
the profiles are smooth and quite well defined. Pressure and 
temperature profiles at this station are shown in Figure 46. 



FIGURE 46 

Pressure and Temperature Profiles at x/rj = 12.4 


The barrel shock is well defined but somewhat smeared. Note the 
smooth monotonic variation of the pressure profile across the 
shear layer and the temperature peak due to viscous (turbulent) 
dissipation. Note also that the pressure monotonically increases 
in traversing the bow shock layer while the temperature remains 
essentially constant due to the effects of bow shock heating (i.e., 
entropy produced) in the curved nose region of the plume. Profiles 
of axial velocity and radial velocity at x/rj = 12.4 are given in 
Figure 47. Note the abrupt change in slope of the radial velocity 
at the viscous dividing streamline location (analogous to that 
obtained in the simple shear layer runs) indicating the diversion 
of streamlines about the hot central core of the mixing layer. 
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FIGURE 47 

Axial and Radial Velocity Profiles at x/rj = 12.4 



5.4.2 Mach 3 Jet into Mach 3 External Stream; Axis to Mixing Layer 
Edge Mapping 


This same calculation was repeated for several inviscid cells using 
a one-domain mapping extending from the axis to the outer mixing 
layer boundary. The inviscid supersonic external flow was repre- 
sented by the use of pressure-flow deflection relations (see 
reference 19) applied at the outer boundary. The calculation was 
performed by SCIPVIS in both the fully-coupled mode and in the 
inviscid limit. The predicted multiple-cell fully-coupled flow 
structure is exhibited in Figure 48 . A comparison of the viscous 
dividing streamline and inviscid slipstream shows them to be in very 
close proximity in the first inviscid cell and then to diverge 
beyond that point with the viscous dividing streamline, remaining 
above the inviscid slipstream. This would indicate that interactive 
effects are fairly weak in the first cell and then increase as the 
width of the mixing region becomes substantial. The axial pressure 
variations shown in Figure 49 confirm this viewpoint. The inviscid 
and fully-coupled Mach disc locations are essentially identical with 
the fully-coupled solution having somewhat greater shock strengths 
due to the additional compression waves generated by the displacement 
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FIGURE 48 

Multiple-Cell Flow Pattern for Underexpanded Mach 3 Jet into Mach 
3 Stream; Single Mapped Domain from Axis to Outer Mixing Layer 
Boundary . 



FIGURE 49 

Inviscid and Fully-Coupled Pressure Variations Along Centerline 
and Outer Boundary. 
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of streamlines about the hot central portion of the mixing layer. 
The inviscid solution exhibits a linear, nondecaying, wave pattern 
beyond the first cell while the fully-coupled solution exhibits a 
nonlinear, mildly decaying wave pattern. 


Profiles of pressure (Figure 50) and temperature (Figure 51) at 
several stations in the first inviscid cell show the gradual build- 
up of viscous interactive effects. At x/rj = 10.6, the two pres- 
sure distributions are nearly identical and the shear layer is quite 
thin. The temperature peaking due to viscous dissipation is not 
pronounced. At x/rj - 20.5, the reflected shock is about to enter 
the plume mixing layer. The two reflected shock positions are 
nearly identical, but the fully-coupled solution has a stronger 
shock. At x/rj = 24.6, the reflected shocks have reached the posi- 
tion of the viscous dividing slipstream or inviscid interface. This 
shock/shear layer interaction is a strongly interactive phenomenon 
and beyond this point, the two solutions start to noticeably diverge 
( see Figure SI). Up until this point , the use of overlaid tech- 
nology 3 is quite reasonable for conditions such as these. However , 
for situations with stronger nearfield interactions as induced by 
faster mixing or chemistry , the overlaid approximation poorly repre- 
sents the plume nearfield. 


5.4.3 Mach 5 Jet into Mach 2 External Stream; Axis to Mixing Layer 
Edge Mapping 


To demonstrate this latter point, the calculation was modified by 
increasing the jet velocity and decreasing the external stream 
velocity keeping all other conditions the same. Here, the jet Mach 
number is boosted to 5 and the external Mach number is decreased to 
2. The nominal velocity ratio has been increased from about 2f/l 
to about 7/1. The predicted flow structure is shown in Figure 52. 
The nearfield interactions are shown to be quite pronounced with the 
viscous dividing streamline and inviscid slipstream positions de- 
parting just beyond the nose region of the plume. The nearfield 
positive displacement effect significantly strengthens the barrel 
shock resulting in a fully-coupled Mach disc position well upstream 
of the inviscid disc position as indicated by the axis pressure 
variation curves in Figure 53 (the fully-coupled disc position is 
at x/rj - 27 while the inviscid disc position is at x/rj ~ 35). 

The subsequent fully-coupled shock locations and cell lengths are 
considerably foreshortened in comparison to the inviscid solution 
due to the substantial curvature of the shocks in traversing the 
highly rotational mixing layer. The shock structure here is signi- 
ficantly stronger than that of the previous Mach 3/Mach 3 case 
(note that the pressures are now plotted on a log scale) since the 
impinging shock waves are largely reflected from the shear layer 
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FIGURE 51 


Comparison of Inviscid and Fully- 
Coupled Pressure Profiles at 
Several Stations. 


Comparison of Inviscid and Fully 
Coupled Temperature Profiles at 
Several Stations. 





















due to the higher vorticity. In the previous case, the ratio of 
reflected to transmitted wave strengths is much less.* 


Profiles at x/r . = 20 shown in Figure 54 exhibit the magni- 
tude of the neaifield interactive effects. The shock propogation 
patterns as evidenced by the pressure and temperature profiles are 
quite different. The fully-coupled barrel shock propogates much 
faster and is at r/rj ~ 1.5 while the inviscid shock is at r/rj ~ 2. 
The cause of the significant displacement effect is quite evident: 
the viscous temperature peak is about 1500 °K above the inviscid 
temperature. This is caused by the very high shear and very high 
Mach numbers in this case. The same interactive effects would 
occur at more realistic velocity ratios and Mach numbers in cases 
with significant afterburning. 


5.5 UNDEREXPANDED JETS INTO STILL AIR 


5.5.1 Mildly Underexpanded , Cold Mach 2 Jet; Pj/Pe = 1* 4 5 


This calculation simulates the experiment of Seiner and Norum 11 * 12 
for a cold (T = 164 °K) Mach 2 jet into still air at a static pres- 
sure ratio (Pj/Pe) of 1.45/1. The predicted flow structure for 
this case (using the kW turbulence model option) is exhibited in 
Figure 55 for the first 40 jet radii in which 8 shock cells occur. 

The nominal upper and lower mixing layer boundaries (where <t> = .05 
and <f> = .95, respectively), the jet mixing layer sonic line, and 
the viscous dividing stream (where ^ = iJ>j) are all depicted. Note 
that the mixing layer engulfs the entire jet downstream of x/rj ~ 

35. The sonic line initially sits above the viscous dividing stream- 
line, intersects it at x/rj ~ 45 and ultimately reaches the axis. 

The details of the predicted sonic line variation for the first four 
cells (Figure 56) are shown to agree quite well with the measure- 


* In the limit of infinite velocity ratio, all the impinging waves 
are reflected. Thus, for a plume into still air, significant 
multiple-cell shock structure exists while for most plumes in 
supersonic flight, shock strengths are usually quite weak beyond 
the first cell. 
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FIGURE 54 

Comparison of Inviscid and Fully-Coupled Pressure, Temperature 
and Axial Velocity Profiles at x/rj = 20. 
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FIGURE 55 

Predicted Flow Structure for Mildly Underexpanded, Cold 
Mach 2 Jet (Pj/P E = 1.45). 
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FIGURE 56 

Comparison of Predicted and Measured Jet Sonic Line Variation 
for First Several Shock Cells. 


82 - 126 


Comparisons of the predicted and measured streamwise pressure vari- 
ations along the jet centerline and along r/rj = .5 are exhibited 
in Figure 57. The comparison is seen to be quite good, indicating 


MACH 2 JET 
PRESSURE VARIATION 



PRESSURE VARIATION 



FIGURE 57 

Comparison of Predicted and Measured Streamwise Static Pressure 
Variations Along Jet Axis and Radial Position r/rj = .5. 
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that the growth of the mixing region (which controls the variation 
of the sonic line and hence the shock cell spacing), and, the 
attenuation of wave intensities by turbulent dissipation, are both 
properly modeled. A comparison of predicted (using kW model) and 
measured u ' 12 longitudinal turbulent intensities is given in Figure 58 
along the jet lipline (r/rj *= 1). Both the measured and predicted 
variations indicate that the turbulence levels peak at the end of 
compression zones and are minimum at the end of expansion zones. 

The comparison is qu ite encouraging in t hat the levels predicted 
based on peak (u^u -2k) and isotroptic (u'u'*=2k/3) estimates of the 
longitudinal intensity bound the measurements (see the discussion 
below) . 



FIGURE 58 

Streamwise Variation of Longitudinal Turbulent Intensity 
Along Jet Lipline. 


&5 



From the discussion in reference 15, we note that the experimental 
data was obtained using a wedge hot-film probe, whose behavioral 
response has been found to differ from that of hot-wires in the 
transonic and supersonic flow regimes. A complete analysis of the 
probe's operational characteristics has been recently published in 
reference 68. However, since only one temperature overheat was 
used to acquire the data for this trial case, the approximate data 
reduction procedure of reference 69 has been used to interpret the 
experimental data presented in this section. 


Since the SCIPVIS code contains a transport equation for the tur- 
bulent kinetic energy, k, an es tim ate to obtain the longitudinal 
turbulence intensity component (u' 2 ) 1/2 /U from k is required. The 
turbulent kinetic energy is defined as: 


k = u’ 2 + v' 2 + w' 2 


From this relationship the longitudinal turbulence intensity can be 
written as: 


(u' 2 )*/U = (2g(x)k)*U -1 


where the function g(x) represents the fractional part of k due to 
longitudinal fluctuations. With shock containing supersonic jet 
plumes, virtually nothing is known about how this function varies 
in the shear layer. Its numerical value can of course never 
exceed 1, and probably never is less than |, which describes a sjtate 
of isotropic turbulence. Actual values taken by the function g(x) 
should fall between these limits, as the predictions indicate (Fig. 


In performing these calculations, a nominal (constant) streamwise 
edge velocity was assigned and the external pressure was set equal 
to the ambient value at all grid points above the sonic line. This 
represents a first approximation to the overall interactive solu- 
tion which adequately represents the flow structure in the first 
several shock cells. Beyond this region, the influence of the ex- 
ternal stream velocity and pressure variations induced by the 
accelerative effects of jet entrainment influence the flow structure 
(see the sensitivity studies of Seiner, et al . 15 exhibiting the in- 
fluence of the small external stream velocity). To account for 


58) . 


96 



these effects, direct coupling with a potential flow solver 16 ' 17 ' 70 
is required. Here, the external flow pressure and velocity pertur- 
bations induced by jet entrainment are obtained by solving the 
potential flow about a geometric body whose surface corresponds to 
the "mean" outer mixing layer edge and along which the source dis- 
tribution is prescribed in accordance with the predicted entrain- 
ment (radial) velocity distribution. The effective geometric body 
and predicted entrainment velocity distribution for this case are 
exhibited in Figure 59. The resultant pressure and streamwise 
velocity variations obtained from the potential flow solution are 
then imposed as outer edge conditions for a revised jet calculation 
and several jet/potential flow iterations are performed until a 
converged solution is arrived at. 



x/rj 


FIGURE 59 

Effective Geometric Matching Surface (teff BODY) for Direct 
Jet/Potential Flow Coupling and Variation of Entrainment 
Velocity (Vent) along Matching Surface. 
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5.5.2 Sensitivity to Turbulence Modeling 


To exhibit the sensitivity of the jet wave field to the representa- 
tion of the turbulent mixing, the calculation above was repeated 
using the ke 31 and ke,cc 39 turbulence models, and was also repeated 
in the inviscid limit. The inviscid prediction (Figure 60) agrees 
reasonably well with the pressure data for the first two shock cells 

MACH 2 JET 
PRESSURE VARIATION 



x/r j 

FIGURE 60 

Comparison of Predicted Streamwise Pressure 
Variations with Data; Inviscid Solution. 
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The damping of wave amplitudes and cell foreshortening effects due 
to the turbulence do not occur in the inviscid limit, and, after 
several shock cells, a linear repetitive pattern persists (subject 
only to numerical damping whose level is negligible in comparison 
to the physical turbulence levels of the problem). The predictions 
made using the "incompressible" ke turbulence model (Figure 61) 
overestimates the rate of turbulent mixing and hence, foreshortens 
the shock cells and dampens the wave intensities prematurely. 

MACH 2 JET 
PRESSURE VARIATION 



FIGURE 61 

Comparison of Predicted Streamwise Pressure Variations with 
Data; ke Solution. 
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The ke,cc prediction (Figure 62) agrees reasonably well with the 
data and appears to slightly underestimate the turbulent mixing 
rate. The performance of the turbulence models in predicting the 
wave structure in this underexpanded jet is in direct accord with 
the centerline decay predictions for the perfectly expanded Mach 
2.2 jet discussed earlier (see Section 5.31 and Figure 37), viz., 
the kW model shows the best overall agreement, the ke model mixes 
too rapidly, and the ke,cc model mixes somewhat too slowly. 

MACH 2 JET 
PRESSURE VARIATION 
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FIGURE 62 

Comparison of Predicted Streamwise Pressure Variations with 
Data; ke,cc Solution. 
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5.5.3 Effects of Jet Temperature and Flight Velocity on Jet Mixing 
and Wave Structure 


The effects of jet temperature and high speed flight velocities on 
jet mixing and wave structure were assessed for the Mach 2 jet into 
still air at the same exit static pressure ratio of 1.45. To assess 
the influence of jet temperature, the jet total temperature was 
boosted to yield an exit temperature of 1500 °K. The predicted flow 
structure for this case is exhibited in Figure 63. The mixing here 


r/r J 



FIGURE 63 

Predicted Flow Structure for Mildly Underexpanded 
(Pj/Pe-1.45) hot (Tj=1500°K) Mach 2 Jet into Still 
Air . 


(using the kW turbulence model) is significantly faster than that in 
the cold jet case with the mixing layer engulfing the entire jet for 
x/rj >25. A comparison of the centerline decay of the species para- 
meter (which for air/air mixing simply represents a tracer species 
and is the same as the velocity decay for a Prandtl number of unity) 
for the two cases (Figure 64) clearly exhibits the differences in 
the mixing rate The streamwise variation in peak turbulent inten- 
sity levels (k 1/2 /U) for the two cases is exhibited in Figure 65. 
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FIGURE 64 

Comparison of <f> Centerline Decay Rates for Cold and 
Hot Jets into Still Air and Hot Jet into Mach 2 Ex- 
ternal Stream. 
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FIGURE 65 

Comparison of Centerline Turbulent Intensities (k ^/U) 
for Cold and Hot Jets into Still Air and Hot Jet into 
Mach 2 External Stream, 


T 





The hot jet has turbulent levels approximately 30 to 40 percent 
larger than those of the cold jet. The predicted variation of cen- 
terline static pressure is exhibited in Figure 66 and indicates a 
more rapid decay of peak pressure levels than in the cold jet case. 


P,atm 



FIGURE 66 

Predicted Centerline Static Pressure Variation for 
Hot, Mildly Underexpanded Jet into Still Air. 


To assess the influence of flight velocity, the same hot jet case 
was run exhausting into a Mach 2 external stream. The predicted 
flow structure is exhibited in Figure 67 where the rate of mixing 
is seen to be considerably slower than the previous two cases due 
to the velocity ratio effect (in the previous two cases, the velo- 
city ratio Uj/U E was essentially infinite; here it is approximately 
2.5/1). The significantly slower mixing rates are also evidenced 
in the species centerline decay (Figure 64), and, the peak turbulent 
intensity levels (Figure 65 ) are seen to be about half of those of 
the hot jet into still air. Of particular interest here is the 
absence of any noticeable wave structure beyond the first shock cell. 
In supersonic/supersonic shear layer interactions (see Section 5.1) 
an impinging wave is partially transmitted and reflected in accord- 
ance with the structural features of the shear layer (whereas for a 
supersonic jet into still air, no transmission occurs and the wave 
reflects from the sonic line). Under the conditions of this case, 
the strength of the reflected wave resulting from the interaction at 
the end of the first shock cell is quite small, and thus, the trans- 
mitted wave intensity is comparable to that of the impinging wave. 
This is evidenced in the predicted streamwise pressure variations 
along the jet axis and along r/rj = .5 (Figure 68) which show no 
noticeable pressure variations beyond the first shock cell. 
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FIGURE 67 

Predicted Flow Structure for Hot, Mildly Underexpanded 
(Pj/Pg=1.45) Jet into Mach 2 External Stream (Uj/Ug~2 .5) 


r/rj=.5 


FIGURE 68 

Centerline and Off-Axis (r/rj=.5) Static Pressure Variation 
for Hot, Mildly Under expanded Jet into Mach 2 External Stream. 





5.5.4 Underexpanded, Cold Mach 2 Jet; Pj/Pe = 3 


To exhibit the treatment of small embedded subsonic regions behind 
Mach discs, the cold Mach 2 jet into still air calculation discussed 
above was repeated with the exit static pressure ratio boosted to 
3/1. The previous case had no Mach discs (regular reflection of the 
shocks at the axis occurred). The predicted flow structure for three 
shock cells using the kW turbulence model is exhibited in Figure 69. 



Predicted Flow Structure for Underexpanded, Cold Mach 2 
Jet into Still Air, Pj/Pe = 3. 


A highly expanded radial scale is utilized to show some details of 
the embedded subsonic domains behind the Mach discs. These domains 
are seen to be quite small for the small radius Mach discs in this 
case. The variation of pressure, Mach number, turbulent kinetic i /2 
energy and scaled turbulent intensity (i.e., u'/U, where u'~(2k/3) 
along the jet centerline are exhibited in Figures 70 and 71. The 
subsonic region is seen to extend about 5 jet radii beyond each disc 
(which corresponds to about 15 Mach disc radii). The mixing process 
is quite rapid and remains confined to the region of the jet center- 
line. The scaled turbulent intensities increase abruptly behind the 
2nd and 3rd discs due to the abrupt decrease in the streamwise velo- 
city in crossing the Mach disc. The turbulent kinetic energy peaks 
in the vicinity of the 2nd Mach disc and then rapidly diminishes, 
peaking slightly in the vicinity of the 3rd Mach disc position. 


105 
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5.5.5 Underexpanded, Hot Mach 2 Jet; Pj/P^ “ 

The hot jet calculations discussed above were repeated with the 
pressure ratio boosted to 5/1 to exhibit other interesting features 
of the flow structure. The predicted nearfield flow structure is 
exhibited in Figure. 72 and a Mach disc of significant radius is 
seen to occur. The jet and Mach disc mixing layers are shaded and 
are shown to intersect at x/rj ~ 15. Note that beyond the Mach disc 



Predicted Nearfield Flow Structure for Underexpanded 
(Pj/P e =5) Hot (Tj=1500 °K) Mach 2 Jet into Still Air. 


position, there are two sonic lines in the flow and for x/rj - 10, 
the code is utilizing all its integration options (i.e., the par- 
tially parabolic option below the Mach disc sonic line; the hyper- 
bolic/parabolic option between the Mach disc sonic line and the 
upper edge of the Mach disc mixing layer, and, between the lower 
edge of the jet mixing layer and jet sonic line; the hyperbolic 
option between the two mixing layers; and, the parabolic option 
above the jet mixing layer sonic line). The overall flow structure 
for 0 < x/rj < 50 is depicted in Figure 73 and encompasses about 5 
shock cells. The mixing engulfs the entire plume at x/rj ~ 35-40. 
The variation of static pressure along the jet centerline is exhi- 
bited in Figure 74. 
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With a Mach 2 supersonic flight velocity, the features of the flow 
structure (Figure 75) change considerably. In particular, no Mach 
disc occurs, the mixing is considerably slower, and the wave struc- 
ture is negligible beyond the second shock cell. This is exhibited 
in Figure 76 where the pressure is seen to be at ambient levels for 
x/rj ~ 30. 
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FIGURE 75 


Predicted Flow Structure for Hot, Underexpanded Mach 2 
Jet into Mach 2 External Stream. 
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FIGURE 76 


Predicted Centerline Static Pressure Variation for Hot, 
Underexpanded Mach 2 Jet into Mach 2 External Stream. 
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5.5.6 Overexpanded, Cold Mach 2 Jet; Pj/Pg = .75 


A typical example of the SCIPVIS code's predictive capability of 
overexpanded flow is shown in Figures 77 and 78 for the Mach 2 
cold jet with P e /P a = 0.749. Figure 77 shows the agreement obtain- 
ed along the jet centerline, and Figure 78 shows that along the 
line r = 0.75 rj. The centerline distribution in Figure 77 repre- 
sents the axial extent of available data for this case, and as is 
evident, the code predicts the irregular looking pressure data with 
remarkable precision. Along the outer radial line in Figure 78, 
there is a good correspondence obtained for the location of the 
first three shock cells; however, the amplitude for the first pres- 
sure maxima differ considerably between measured and predicted 
results. It is known that the oblique shock strength is a strong 
function of radial position as the sonic line is approached, so 
that small differences in plume spreading could account for this 
discrepancy in the amplitude of the pressure maxima. Figure 78 
also indicates that there is a marked difference in the shock wave- 
length of the downstream shock structure, again suggesting that 
the measured jet plume has a faster plume spreading rate. 



FIGURE 77 

Comparison of Predicted (kW) and Measured (References 11 & 12) 
Streamwise Pressure Variation along Jet Centerline for Mach 2 
Jet . 
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FIGURE 78 

Comparison of Predicted (kW) and Measured (References 11 & 12) 
Streamwise Pressure Variation along R/Rj = .75 for Mach 2 Over- 
expanded Jet. 


5.5.7 Underexpanded, Cold Mach 1.4 Jet; Pj/Pg = * • 36 


For this exit Mach number, only one trial case is considered; the 
mildly underexpanded flow where Pe/Pa = 1*36. In these comparisons 
the experiment data 11 ’ 12 is shifted downstream .233 rj to correspond 
to a position where the static ports are located on the probe body. 

With lower Mach number flows the static probe port position provides 
a more realistic account of stream Mach number when used in conjunction 
with a supersonic total pressure tube. This was not required for the 
Mach 2 condition, since the static probe is designed for this elevated 
Mach number. 


Figures 79, 80, and 81 respectively show the comparison between the 
measured and predicted static plume pressure for the ke2 , ke2,cc 
and kW turbulence closure models. The upper part of each figure 
displays the centerline results, which represents the axial extent 
of measured data, and the lower part of each figure the results 
along the radial line r = 0.9 rj. It is quite evident from all 
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FIGURE 79 

Comparison of Predicted (ke2) and Measured Streamwise Pressure 
Variations for Mach 1.4 Underexpanded Jet. 
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centerline comparisons for this case that all three turbulence clo- 
sure models predict the amplitude of the measured data extremely 
well, but that the ke2— cc model slightly overpredicts the location 
of the pressure maxima. This can be more readily observed in the 
off-axis data which clearly shows that the ke2-cc model overpredicts 
and the ke2 model underpredicts the wavelengths of the downstream 
shocks. As was found above for the Mach 2 case, the kW model ap- 
pears to provide the best match to experimental observations. How- 
ever, it is not as yet entirely clear that the kW option would work 
best when the potential flow coupling is incorporated to account for 
entrainment effects on the external stream. 


5.5.8 Underexpanded, Cold Sonic Jet; Pj/Pg = 1«62 


Similar to what was done in the Mach = 1.4 case, the experimental 
data has been shifted in these comparisons 0.241 rj downstream to 
correspond to the locations of the probe's static ports (the sonic 
nozzle has a smaller exit diameter). The results shown in this sec- 
tion were made for the underexpanded condition Pe/Pa = 1.62, which 
matches previously acquired experimental data 11 ’ 12 for the sonic noz- 
zle. Figures 82 and 83 demonstrate respectively the agreement 
obtained between the measured and predicted static plume pressure 
for the ke2 and ke2-cc turbulence closure models. The upper part 
of each figure is for the centerline data, whereas the lower part 
of each figure is for the off-axis measurements and predictions. 

The measurements in the off-axis data refer to the radial location 
r = 0.9 rj, and the predictions to the radial location r = 0.75 rj. 

By virtue of the agreement between the first few pressure maxima, 
this difference in radial location suggests that there is a signifi- 
cant difference between the predicted and measured plume structure. 
Also unusual is the very rapid decay of the measured static pressure 
distribution along the centerline after the fifth shock cell. This 
uncharacteristic behavior strongly suggests that new experimental 
information be acquired to examine this large difference in pre- 
dicted and measured results. At this time, it is not known to what 
extent the predicted results depend on the selection of a "nominal" 
first-pass external flow Mach number used to simulate ambient con- 
ditions. Comparing Figures 82 and 83 shows that the ke2 model 
matches the measured data better than the ke2-cc model, for the con- 
ditions chosen. (The kW turbulence model results for this case 
fall between the ke and ke-cc results and are not exhibited). This 
sensitivity is further discussed in reference 15 and will be resolved 
when potential flow coupling with SCIPVIS is available to provide 
"proper" external flow conditions. 

















6.0 


CONCLUDING REMARKS 


A new two-dimensional, interactive jet mixing model, SCIPVIS, 
has been described which provides for the fully-coupled analysis 
of under or overexpanded jets exhausting into supersonic or 
quiescent external streams. SCIPVIS combines shock-capturing 
methodology for analyzing supersonic mixing regions with simple 
pressure-split methodology for analyzing subsonic mixing regions. 

A direct coupling approach has been introduced for the interactive 
coupling of the jet and external flow solutions along the 
outer jet boundary. This coupling approach has entailed mapping 
the jet boundary onto a constant mapped coordinate surface and 
developing viscous/inviscid growth rules for the jet boundary 
variation. This coupling approach is directly extendable to 3D 
jets and is ideally suited for coupling with a subsonic/transonic 
potential flow or Euler solver. This iterative coupling entails: 

(1) defining a matching surface which encloses the viscous/inviscid 
jet; (2) solving the jet flowfield with the pressure and streamwise 
velocity variation imposed on this surface as outer boundary 
conditions; and (3) solving the external flowfield with an inner 
transpiration (source/sink) boundary condition imposed on this 
surface corresponding to the entrainment (normal) velocity 
distribution determined from the jet solution. 

The ability of SCIPVIS to treat a variety of ducted and free jet 
mixing problems has been demonstrated via the performance of 
numerous calculations. SCIPVIS can readily handle the full range 
of under or overexpanded conditions encounterable in conventional 
aircraft exhaust problems. Its application to highly underexpanded 
flow problems as encounterable in the analysis of rocket and 
ramjet exhausts requires that the parabolization of the flow equations 
in the highly inclined nearfield shear layer region be performed 
in a plume-slipstream-oriented coordinate system. This parabolization 
is most readily accomplished via the use of a spherical coordinate 
system with the origin selection based on the initial plume slip- 
stream angle. 

Comparisons of predictions with laboratory data for mildly under 
or overexpanded jets exhausting into still air have indicated that 
SCIPVIS has the capability of accurately reproducing the detailed 
multiple-cell wave structure observed including the attenuation 
of wave intensities and curvature of wave fronts due to wave/ 
turbulent interactions. This unique capability can pave the way 
for extending SCIPVIS into a program for estimating shock noise 
in supersonic jets 15 . Immediate problem areas to be resolved 
include the influence of entrainment effects on the "farfield" 

(last several shock cells) wave structure 15 for a quiescent 
environment, and, the selection of an appropriate turbulence model. 
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The coupling of SCIPVIS with a potential flow code will lead to 
resolution of this first problem area. Our present analysis of 
the CALSPAN supersonic jet mixing data 65 should lead to 
improvements in two-equation turbulence models which permit them 
to reliably analyze the range of jet flow conditions encounterable 
in aircraft exhausts. 

The analysis of 2D or 3D jet/af terbody interaction problems in 
a subsonic/transonic external flow environment requires the 
development of sophisticated pressure-split or velocity-split 
methodology to treat the subsonic portion of the jet mixing layer. 
In the simplest of rectangular jet mixing problems, cross flow 
vortices can be generated in the corner regions whose prediction 
requires a robust elliptic cross flow integration procedure. 
Whereas the shock-capturing methodology in SCIPVIS can readily 
be extended to treat complex supersonic 3D jet mixing processes, 
the simple pressure-split methodology in SCIPVIS utilized in 
subsonic zones, cannot. To address the generalized subsonic jet 
mixing problem, a separate model, SPLITP, has been developed which 
utilizes the same mapped computational network, boundary growth, 
and coupling procedures as SCIPVIS. However, SPLITP integrates 
the nonconservative form of the parabolized jet mixing equations 
utilizing an implicit, upwind algorithm. In subsonic regions, 
the cross flow velocity and pressure variations are determined 
using elliptic methodology. Weak wave processes in supersonic 
regions are analyzed by a local, explicit characteristic type 
solver now under development. Thus, SPLITP comprises a parabolic/ 
elliptic model with limited hyperbolic capabilities while SCIPVIS 
comprises a parabolic/hyperbolic model with limited elliptic 
capabilities. The unification of these two models into a single 
code will be an ultimate goal of this program. Since both models 
employ the same variables, mappings, turbulence models, etc., this 
unification process is straightforward, viz. only the streamwise 
and cross-flow integration subroutines in the two codes differ. 
SPLITP has already been extended to analyze curved boundary layer 
flows 16 ’ 17 and thus, the combined SCIPVIS/SPLITP model can be 
used to analyze the (nonseparated ) afterbody boundary layer and 
nozzle flow regions as well as the plume flowfield. A complete 
description of SPLITP and the steps entailed in integrating SPLITP 
and SCIPVIS into a unified model will be given in Part II of this 
report . 

Significant progress has already been made in extending SCIPVIS 
to analyze 3D exhaust problems. A preliminary version of this 
model is operational for supersonic rectangular nozzle/exhaust 
problems and many of the earlier inviscid results of Dash and 
DelGuidice 7 1 ’ 7 2 have been reproduced and/or improved upon. A 
series of fundamental 3D wave/shear layer interaction numerical 
studies are in progress akin to the 2D studies described in this 
report, and various mappings/coordinate networks are being 
investigated to treat a range of exhaust nozzle configurations. 
This work will be documented in the near future. 
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